A perspective on the interfacial properties of nanoscopic liquid drops 
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The structural and interfacial properties of a nanoscopic liquid drops are assessed by means of 
mechanical, thermodynamical, and statistical mechanical approaches that are discussed in detail, 
including original developments at both the macroscopic level and the microscopic level of density 
functional theory (DFT). With a novel analysis we show that purely macroscopic (static) mechanical 
arguments can lead to a qualitatively reasonable description of the surface tension and the Tolman 
length of a liquid drop; the latter parameter which characterizes the curvature dependence of the 
tension is found to be negative and has a magnitude of about a half of the molecular dimension. 
A mechanical slant cannot, however, be considered satisfactory for small finite-size systems where 
fluctuation effects are strong. From the opposite perspective, a curvature expansion of the macro- 
scopic thermodynamic properties (density and chemical potential) is then used to demonstrate that 
a purely thermodynamic approach of this type can not in itself correctly account for the curvature 
correction of the surface tension of liquid drops. We emphasize that any approach, such as, e.g., 
classical nucleation theory, which is based on a purely macroscopic viewpoint does not lead to a 
reliable representation when the radius of the drop becomes microscopic. The description of the 
enhanced inhomogeneity exhibited by small drops (particularly in the dense interior) necessitates a 
treatment at the molecular level to account for finite-size and surface effects correctly. The so-called 
mechanical route which corresponds to a molecular-level extension of the macroscopic theory of elas- 
ticity, and is particularly popular in molecular dynamics simulation, also appears to be unreliable 
due to the inherent ambiguity in the definition of the microscopic pressure tensor, an observation 
which has been known for decades but is frequently ignored. The union of the theory of capillarity 
(developed in the nineteenth century by Gibbs and then promoted by Tolman) with a microscopic 
DFT treatment allows for a direct and unambiguous description of the interfacial properties of 
drops of arbitrary size; DFT provides all of the bulk and surface characteristics of the system that 
are required to uniquely define its thermodynamic properties. In this vein, we propose a non-local 
mean-field DFT for Lennard- Jones (LJ) fluids to examine drops of varying size. A comparison of 
the predictions of our DFT with the recent simulation data based on a second-order fluctuation 
analysis [J. G. Sampayo et ah, J. Chem. Phys. 132, 141101 (2010)] reveals the consistency of the 
two treatments. This observation points out the significance of fluctuation effects in small drops, 
which give rise to additional entropic (non-mechanical) contributions, in contrast to what one ob- 
serves in the case of planar interfaces which are governed by the laws of mechanical equilibrium. 
A small negative Tolman length (which is found to be a tenth of the molecular diameter) and a 
non-monotonic behaviour of the surface tension with the drop radius are predicted for the LJ fluid. 
Finally, the limits of a validity of the Tolman approach, the effect of the range of the intermolecular 
potential, and the behaviour of bubbles are briefly discussed. 



I. INTRODUCTION 

The study of inhomogeneous systems presents a much 
more significant challenge than that of homogenous flu- 
ids. By definition the non-uniformity of the number den- 
sity throughout the sample adds a mathematical com- 
plexity to the theoretical description - the correlation 
functions become multivariable functions and, within a 
variational formalism, the thermophysical functions be- 
come functionals of the single-particle density, and as 
a result the partial derivatives relating the equilibrium 



properties with a particle density must be replaced with 
the corresponding variational expressions. At the level of 
a formal physical description, a more fundamental issue 
arises: the thermodynamic quantities that are familiar 
in studies of uniform fluids can not always be defined 
properly in the inhomogeneous region. The interface be- 
tween vapour and liquid phases or two liquid phases are 
ubiquitous examples of non-uniform systems. The study 
of two bulk phases separated by a planar interface (and 
stabilized by an arbitrary weak external field) does not 
present a particular problem, since the non-uniqueness 
in the definition of the local variables (such as the pres- 
sure tensor or the position of the interface) does not give 
rise to an ambiguity in the measurable quantities (such 
as the surface tension). The situation is not. however. 
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as straightforward for systems exhibiting spherical sym- 
metry, such as small liquid drops (liquid surrounded by 
vapour) or bubbles (vapour surrounded by liquid), and 
our goal in this paper is to provide a critical discussion 
of the most popular methodologies for treating the inter- 
facial properties of such systems. 

There are three general routes to the determination 
of interfacial properties of a small fluid droplet or bub- 
ble: these involve the choice of a mechanical, a thermo- 
dynamical, or a statistical mechanical description. The 
first successful mechanical description dates back to the 
beginning of nineteenth century when Young [T] and 
Laplace [2] derived a relationship for the difference in 
pressure p between a phase a on one side of a curved in- 
terface and the surrounding phase f3. For a macroscopic 
system their simple relation can be expressed as 



Pa - Pp 
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where 7 is the surface tension playing the role of the 
restoring force acting against changes in the area of the 
interface, and R is the radius of the drop or bubble; 
though Young's paper predates Laplace's more thorough 
derivation by a few months, the expression is more com- 
monly referred to as the Laplace equation probably be- 
cause Young only describes the dependence of curvature 
in words and not as an explicit formula [3J . Both deriva- 
tions rely on a macroscopic definitions of the quantities 
p a , pp, and R and as that they are all assumed to be 
uniquely defined. In particular, one assumes that both 
phases behave in the same way as the corresponding bulk 
phase. However, when one considers smaller and smaller 
drops, surface contributions propagate progressively into 
the interior of the drop so that the density profile be- 
comes highly structured (as we will show later in our 
paper), and the concepts of the "bulk" density of the liq- 
uid, the scalar pressure, and the radius of the drop lose 
their obvious characteristics; for an in-depth review of the 
problematic issues associated with spherical surfaces see 
the excellent account by Henderson [4] . A natural exten- 
sion of the concept of the scalar pressure to non-uniform 
systems is an introduction of a second-rank tensor P(r) 
a local quantity related to the force between the inter- 
acting molecules at a point r. In the absence of external 
fields, the sum of all forces P(r) • h&A acting on the in- 
finitesimal area AA, where n is the unit vector normal to 
the particular element of area, must balance: 



/ P(r)-racL4 
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As is customary, Gauss's divergence theorem can be 
employed to re-express the surface integral as one in 
the divergence of the pressure tensor over volume V. 
The resulting equilibrium condition must apply to each 
infinitesimal element so that a microscopic mechanical 
treatment then relies on the simple condition for mechan- 
ical equilibrium at every point in the system [3J: 



It is clear, however, that Eq. ^ cannot be used to define 
the pressure tensor uniquely because any tensor P'(v) 
which differs from P(r) by a curl still satisfies the equi- 
librium condition (J3J) . Even though the surface tension, 
which, for a planar interface, can be obtained from the 
pressure tensor as 



7= [P n (z) - P t (z)]dz : 



(4) 



where P n {z) = Pbuik and Pt(z) are the normal and tan- 
gential components of P(r), is invariant to the particular 
form chosen for the local pressure tensor, the first mo- 
ment of the difference between the two components, 
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z[P n (z) - P t (z)]dz , 



(5) 



is not [3J 03 [5] . The latter defines the so called surface 
of tension, i.e., the surface at which the surface tension 
is deemed to act; the surface of tension plays a crucial 
role in the determination of the curvature dependence of 
the surface tension as will be discussed in the subsequent 
discussion. 

The first thermodynamic expression relating bulk ther- 
modynamic properties with the ones associated with the 
system of a liquid drop (at the same temperature T) is a 
familiar Kelvin equation [7] 
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where p v (R) is the vapour pressure of a drop of radius 
R (i.e., p^, at is the saturation pressure), 700 is the sur- 
face tension associated with the planar interface, and kg 
is the Boltzmann constant. The thermodynamic route 
to spherical interfaces has been further developed by 
Gibbs [HHIO] and Tolman [TTHIl] (with further develop- 
ments by Buff [5], Koenig [14] , Hill [TS], and Kondo [TB] 
amongst others). According to Gibbs [5J the curva- 
ture corrections become essential only for very small 
droplets, with the conjecture that the surface tension 
decreases monotonically on decreasing the radius of the 
droplet. From other perspectives the Thomsons (father 
and son) [17] suggested that one should allow for the 
possibility of a non-monotonic behaviour (actually a min- 
imum in the surface tension with decreasing radius fol- 
lowed by one or more maxima) , while Bakker |18) insisted 
on the invariance of the surface tension of the droplet 
with its radius. Tolman |13j . who can be considered as 
one of the main promoters of the Gibbsian view, put for- 
ward a rigorous theory for the dependence of the surface 
tension on the radius of the drop based purely on ther- 
modynamic arguments (essentially assuming that high- 
frequency, short-wavelength, capillary wave terms and 
elastic deformations from spherical geometry are negli- 
gible as one would expect at the level of leading order in 
curvature [I]): 



V • P(r) =0. 



(3) 
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Here, j(R) is the surface tension of the droplet of an ar- 
bitrary radius R, and 5 = R e — R s is the difference in the 
distance between the surface of tension R s (where the 
tension acquires its minimum) and the Gibbs equimolar 
dividing surface R e (where the excess superficial density 
of particles effectively vanishes); the so called Tolman 
length is that corresponding to the value of S(R) in the 
limit of the planar interface S = limR-^ 5(R) = z e — z a , 
with the appropriate distances from the interfacial plane 
now represented by z e and z s . According to the Gibbs- 
Tolman view of a decrease in the surface tension with 
decreasing radius, the Tolman length would thus be a 
positive quantity, S > 0. It is important to realize that 
while for planar interfaces the choice of dividing surface 
is arbitrary, it is a "necessity, not merely a convenience" 
for systems with curved interfaces [3]. Owing to the 
phenomenological origin of thermodynamic approaches of 
this type, the description is expected to become increas- 
ingly inappropriate when one is attempting to represent 
smaller and smaller droplets as was pointed out early 
on by Farkas [19], by Guggenheim [20], and by Tolman 
himself [13] . Though Tolman also incorporated higher 
order terms in the radius dependence, it is question- 
able to what extent these are meaningful (or even phys- 
ically relevant) for very small droplets; this is because, 
to higher order in curvature, the value of the surface 
tension becomes dependent on the choice of the divid- 
ing surface [311]. Before proceeding we should, however, 
acknowledge that the description of curvature deforma- 
tions for non-spherical geometries beyond the first-order 
correction are in common use, particularly in treating 
complex fluids with low tensions such as surfactant ag- 
gregates and membranes (e.g., see Refs. [2TH24] ) : this in- 
troduces additional complications which are beyond the 
scope of our paper. 

Apart from the aforementioned issues with the me- 
chanical or thermodynamical treatment of curved inter- 
faces, neither of the approaches provide us with a direct 
link between the microscopic (intermolecular interactions 
and local structure) and macroscopic (thermodynamic) 
properties of the fluid. This is possible with a full sta- 
tistical mechanical description of the non-uniform fluid. 
Classical density functional theory (DFT) provides one 
with a very powerful tool for the description of inho- 
mogeneous systems P5T - [2"T] . In a purely mechanical or 
thermodynamical approach one manipulates local many- 
body quantities such as the force, local energy, local pres- 
sure etc. to describe the interfacial properties, but these 
can often be ill-defined as there is no unique way of as- 
signing contributions from the intermolecular forces to a 
particular element of space [25]. By contrast, in a DFT 
treatment the full partition function (and therefore ther- 
modynamic potential) of the system is formulated ex- 
plicitly in a spatially dependent form in terms of the sin- 
glet density which is a well-defined one-body function, 
allowing for a unique description of the thermodynamic 
properties. Such an approach is, however, still not en- 
tirely straightforward for inhomogeneous systems char- 



acterized by curved interfaces [U , and care has to be 
taken with the precise route that one employs to compute 
the interfacial properties. In its original form, the DFT is 
formulated in the grand canonical ensemble in which an 
isolated finite-size drop of liquid is unstable with respect 
to its vapour. This leads one to an inevitable key ques- 
tion: How does one stabilize a drop of fluid of finite size? 
Assuming that a stabilized drop can then be examined to 
determine the equilibrium density profile by minimizing 
the appropriate functional, one then has to establish a 
unique and consistent methodology for the desired ther- 
modynamic and interfacial quantities from a knowledge 
of structure of the fluid. 

The system of an isolated drop of liquid surrounded 
by its vapour (or the inverse case of an isolated bubble 
of vapour in a liquid) is ubiquitous and has been studied 
extensively by experiment, theory and molecular simu- 
lation since the early description of Young and Laplace. 
Experimental evidence of the effect of the system size 
and the curvature on the surface tension is scarce as the 
variation from the macroscopic (planar) values is directly 
measurable for only very small dimensions. Early indi- 
cations of the curvature effects were obtained by Reinold 
and Riicher |30j from experiments on thin films of soap 
solutions, the thickness of which can be estimated from 
the colour: the surface tension was found to be con- 
stant down to a thickness of ~ 50 nm, followed first by 
a decrease with decreasing thickness and then by an in- 
ferred increase. Though this lead to the Thomson view of 
the possibility of a non-monotonic behaviour in the ten- 
sion [17] , such an analysis should be made with particular 
care owing to the inherent difference between aqueous so- 
lutions of amphiphilic compounds (surfactants which will 
accumulate at the interface to differing degrees, and lead 
to a decrease in tension) and fluid drops of pure sub- 
stances. At the turn of the twentieth century Weber [31] 
also detected evidence of size effects on interfacial prop- 
erties in his experiments of the contact angle in oil- water 
systems, though again the findings are difficult to inter- 
pret in such mixtures. In the more recent analysis of 
experimental data for the effect of curvature in fluid sys- 
tems, the Tolman relation is often employed at leading 
order in curvature, together with the measured vapour 
and liquid densities and the surface adsorption, to esti- 
mate the surface tension: for example, in the case of a 
drop of water the surface tension is found to remain es- 
sentially constant (within a few percent) for radii down 
to ~ 10 nm, and to decrease rapidly thereafter [32]). 
However, as already been emphasized, for small drops 
one is at the limit of the applicability of thermodynamic 
approaches of this type, and any tautological conclusions 
of this kind should be viewed with some scepticism. This 
having been said, a macroscopic treatment continues to 
be in use without reservation to this day, e.g., see the 
work of Xue et al. [33] . The surface force apparatus was 
used early in its development by its pioneers to provide a 
direct measure of the interfacial forces of curved surfaces 
at the microscopic level: Fisher and Israelachvili [53H31)] 
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investigated the limit of validity of the Laplace and Tol- 
man relations (using the corresponding Kelvin [7 j macro- 
scopic thermodynamic description of the vapour pressure 
of curved interfaces) for a meniscus of hydrocarbon fluid 
between mica spheres/cylinders; notwithstanding some 
complications due to impurities, the Kelvin relation is 
found to be valid for menisci with radii down to ~ 4 nm 
(corresponding to about 10 diameters of typical small 
molecules) , with a marginal possible improvement in the 
description of the data for a Tolman-like dependence of 
the tension with curvature. We should note however that 
the analysis of Fisher and Israelachvili at very high cur- 
vature (small radii) was brought into question in a later 
study by Christenson [37], and the effect has now been 
found to be very sensitive to differences in the structure 
and polarity of the molecules [55] ■ A more recent lattice 
gas Monte Carlo study of the atomic force microscope 
experiment has also indicted that there is a lower limit 
in the size of the system (corresponding to radii of about 
2 nm) below which it is no longer possible to stabilize a 
meniscus of fluid [39] , Our overall understanding is not 
helped by the analysis of data for the deformation of fluid 
interfaces obtained from small angle X-ray and neutron 
scattering experiments, which for the vapour-liquid inter- 
face of water and organic molecules j40j [41] is consistent 
with a negative Tolman length, S < 0, while in the case 
of surfactant monolayers [12] 25J the data supports the 
original Gibbs- Tolman picture with S > 0; care should 
again be taken with the analysis for the more complex 
systems comprising amphiphilic compounds. 

The body of work on molecular simulation of vapour- 
liquid drops and bubbles though extensive is understand- 
ably not as sizeable as that for its planar counterpart 
(see Ref. [H] for a recent review of the latter). In 
one of the first continuum studies of liquid drops car- 
ried out a few years before the better known work of 
Binder and co-workers [45] 06] with lattice gas models, 
Rusanov and Brodskaya [47j examined drops of trun- 
cated Lennard- Jones (LJ) particles inside a spherical 
hard cavity by molecular dynamics (MD) simulation, 
and calculated the pressure tensor of the system. Ru- 
sanov and Brodskaya showed that one cannot obtain a 
uniform value of the tensorial components of the pres- 
sure in the centre of small drops (bringing into ques- 
tion a validity of the macroscopic mechanical definition, 
though admittedly the uncertainty in the computed val- 
ues is large) , and instead calculated the tension from the 
Laplace relation with the pressure of the liquid interior 
obtained in a thermodynamically consistent way from a 
bulk system with an equivalent chemical potential. In 
agreement with the Gibbs- Tolman view, the tension was 
found to decrease with decreasing drop radius. Powles et 
al. E5H3H] also simulated drops of the LJ fluid (essentially 
for the full range of the pair interaction) in coexistence 
with vapour using standard periodic boundary conditions 
and determined the tension and Tolman length from the 
Kelvin thermodynamic relation. In their well cited pa- 
per, Thompson et al. |50) reported values of the surface 



tension of shifted and truncated LJ drops within soft- 
wall cavities and with dynamic walls (tied to the centre- 
of-mass of the drop) obtained by MD simulation from 
both mechanical (pressure-tensor) and thermodynamic 
(Tolman and Laplace) routes; Thompson et al. acknowl- 
edged the problems associated with the use of such ap- 
proaches to determine the surface tension for small drops, 
and highlighted the inadequacy of the Laplace and Kelvin 
relations for drop radii smaller than about 10 molecular 
diameters. These early simulation studies all appear to 
confirm the Gibbs- Tolman view of a decrease in the sur- 
face tension with decreasing drop size, corresponding to 
a positive value of the Tolman length, 6 > 0; one should 
stress, however, that only relatively small systems were 
examined, and the Tolman length is certainly expected 
to depend on the system size and on the range of the 
interactions (as we show later in paper). 

There have since been a number of computer simula- 
tion studies of liquid drops [5THT9"] and bubbles [5DH59"] ; 
here we refer to some of the representative work where 
new findings relevant to our current study are reported, 
and we make no attempt to provide a full review of all the 
literature on curved interfaces of mixtures, nucleation, 
cavitation, and other non-equilibrium processes. In one 
of the first large-scale simulation studies, Nijmeijer et 
al. [52] indicated that due to statistical scatter there is 
a large uncertainty in the sign of the Tolman length; the 
main finding being that its value is close to zero, a con- 
clusion supported by one of the most recent simulation 
study [79] . The small absolute value of the Tolman length 
is generally supported by the more recent simulation data 
(e.g., Refs. [55] ESI EH [75]) for drops of up to - 8 x 10 5 
LJ particles corresponding to radii of almost 100 diame- 
ters [62] , It is also very clear that the fluid drops experi- 
ence large fluctuations in shape and size, particularly in 
the case of large systems, as has been shown by Arcidia- 
cono et al. [61] and Salonen et al. [64]. With a thermo- 
dynamic approach based on a linear response of the free 
energy to small volume perturbations, El Bardouni et 
al. [55j reported some values of the tension for spherical 
and cylindrical surfaces that are larger than the planar 
limit (corresponding to S < 0), though the uncertainty is 
such that they concluded that the tension is essentially 
curvature independent. This is not the case in the studies 
carried out by Vrabec et al. 63 , who used the conven- 
tional mechanical (pressure-tensor route) and found that 
the surface tension decreased sharply and monotonically 
with decreasing drop radius (5 > 0). In related stud- 
ies of nucleation in fluids ten Wolde and Frenkel |54j . 
and Neimark and Vishnyakov |67j have shown that er- 
roneous nucleation barriers result from the use of pres- 
sure tensors (mechanical route), and that the Tolman 
equation is not valid for clusters with radii below four 
molecular diameters [57]. More recently, Binder and co- 
workers et al. [Ml IM E2 [23 [78] have used a thermo- 
dynamic analysis using a Landau free energy and grand- 
canonical Monte-Carlo approach to determine the sur- 
face free energy and interfacial tension for drops of vary- 
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ing size, finding that the curvature corrections cannot 
be described with the simple Tolman relation for small 
drops. These authors also find that the interfacial ten- 
sion increases above that of the planar interface, albeit 
very marginally, for drops with radii larger than about 
8 molecular diameters, which points to a small and neg- 
ative Tolman length. Though this finding is in contra- 
diction with the large body of work based on a purely 
mechanical analysis of the simulation data, it is consis- 
tent with the earlier study of El Bardouni et al. |55j . 
with a thorough analysis based on the Laplace relation 
for very large drops [71] and with the use of test-area 
deformations [72] . We shall come back to this interesting 
feature later in our discussion. 

In the case of bubbles within an fluid of LJ particles, 
Park et al. [82] have used the mechanical expression for 
the normal and tangential components of the pressure 
tensor and the Laplace relation to estimate the tension 
and Tolman length; they find that though the tension of 
the bubble is now greater than that of the planar inter- 
face, the Tolman expression for the first-order curvature 
correction does not quantitatively reproduce the calcu- 
lated surface tension of the bubble, possibly due to an 
inconsistency in the calculation of the Tolman length. 
By contrast, in their recent study of very small LJ bub- 
bles, Matsumoto and Tanaka [55] determined the vapour 
pressure with an empirical equation of state (rather than 
via a pressure-tensor route) , finding that the surface ten- 
sion is independent of the radius of the bubble (which 
corresponds to S = 0) , and confirming the validity of the 
Laplace relation for radii down to ~ 1.7 nm (in terms of 
the LJ parameters for argon). However, the latest esti- 
mates of the Tolman length for bubbles by Block et al. 
[73"] now suggest a small negative Tolman length (corre- 
sponding to about a tenth of the molecular diameter) as 
in the case of liquid drops. 

In view of the disparate findings reported in the various 
simulation studies of fluid drops and bubbles it would not 
be unfair to say that there is still no clear consensus re- 
garding the curvature dependence of the surface tension 
and the sign of the Tolman length. Different (essentially 
macroscopic) routes are employed to analyze the data 
for the interfacial properties, the validity of which are in 
question for small systems. To add to the confusion the 
treatment of the range of the intermolecular potential 
(long ranged, versus truncated or truncated and shifted 
potentials) has been the bane of the calculation of the 
surface tension for planar interfaces (particularly in ap- 
proaches employing a mechanical route because of the 
discontinuous nature of the forces), leading to general 
conclusions which are in apparent conflict; see the pa- 
per by Trokhymchuk and Alejandre [90] and references 
therein. The contradictory data for curved interfaces are 
most certainly also compounded by the treatment of the 
range of the potential, as Lei et al. [62] have demon- 
strated for large liquid drops, reiterating the fact that 
the surface tension and Tolman length are very sensitive 
to the value of the intermolecular potential cutoff that is 



employed. 

The full armoury of phenomenological thermodynamic 
approaches and the more sophisticated statistical me- 
chanical theories have been employed to describe the 
interfacial properties of systems with curved interfaces, 
including mean-field, square-gradient (generalized van 
der Waals), capillary- wave, density functional, and fun- 
damental measure theories [91-160]. In the following 
discussion we will again not focus on studies of nucle- 
ation or criticality, which represent entire fields in them- 
selves. The general conclusions that can be drawn from 
these studies are as inconclusive as those gleaned from 
direct molecular simulation. As we have already men- 
tioned, using his macroscopic thermodynamic approach 
Tolman [13] found a monotonically decreasing surface 
tension with decreasing drop radius, corresponding to 
5 > which was of the order of 0.1 nm; if one ex- 
tends the concept of the Tolman length to a function 
S(R) = Re — R s of the drop radius then a non-monotonic 
dependence of the surface tension with the radius can be 
obtained [150] . 

Hemingway et al. [108] have compared thermody- 
namic, mechanical, and statistical mechanical routes for 
vapour-liquid surface tension and Tolman length of the 
penetrable-sphere model. This provides evidence of the 
consistency between the thermodynamic and statistical 
mechanical routes (though it cannot be considered as a 
proof), while in the case of a mechanical treatment the 
value of the Tolman length depends on the choice of lo- 
cal pressure tensor (as first demonstrated by Schoficld 
and Henderson [5] and later by Blokhuis and Bedeaux 
[125] 1 . One can formulate a form of the local pressure 
tensor that gives a unique expression for the surface of 
tension in the case of systems with spherical symmetry 
as demonstrated by Baus and Lovett [1201 11241 1126] , but 
the expressions are much more complicated and there are 
issues in their implementation to liquid drops [124] 1128] . 
Sampayo et al. [72] have also shown that a virial re- 
lation only corresponds to the leading-order term in the 
free-energy change due to the deformation of small drops, 
and that there are large contributions from the second- 
order (fluctuation) term with a magnitude which is com- 
parable to the first-order contribution. Lekner and Hen- 
derson [104] have demonstrated that the first-order con- 
tribution to the change in free energy accompanying a 
change in the interfacial area captures the entire mechan- 
ical contribution that one would obtain for the difference 
in the appropriate components of the pressure tensor (cf. 
the Irving-Kirkwood [94] expression in the case of a pla- 
nar interface). It is therefore clear that first-order me- 
chanical routes which rely on the pressure tensors are to 
be avoided for small drops as they do not incorporate 
the large contributions due to thermal fluctuations. The 
main advantage of the penetrable-sphere model is that it 
can be solved exactly at the mean-field level at zero tem- 
perature where Hemingway et al. [108] find a negative 
Tolman length, 5 = —cr/2, with a the molecular diame- 
ter. It is not clear that relation will still hold at higher 
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temperatures, however, the main problem with such an 
approach is the lack of knowledge of a good approxima- 
tion for the direct correlation function for generic fluid 
models [TUB]. 

The square gradient theory (SGT), which belongs to 
a class of more general density functional theories |25j . 
is rooted in van der Waals' [92] original treatment for 
fluid interfaces (and in the earlier work by Rayleigh [9"T]). 
which was rediscovered and popularized by Cahn and 
Hilliard [321 ■ Before we discuss the findings of mi- 
croscopic SGT approaches, we should briefly mention 
the related phenomenological treatment referred to as 
capillary-wave theory. As Henderson 4J has pointed out 
an analogy with hydrodynamics can be made to exam- 
ine the surface tension of a fluid as the restoring force 
due to thermally excited surface waves; frequent use of a 
capillary- wave description has be made to describe planar 
interfaces |100|, I102|, I103J and to represent liquid drops 
[5[ 1 108] . The bare capillary- wave surface tension cor- 
responds to the equilibrium (infinite wave-length free- 
energy contribution) thermodynamic surface tension in 
the case of a planar interface and also to that of a spher- 
ical interface at the level of leading order in curvature 
[H I121j . This means that the surface tension of the 
system is required as an input if one wants to employ 
capillary wave approaches to describe interfacial systems. 
One of the first to use SGT to examine curved interfaces 
and liquid drops were Falls et al. [107j : they approxi- 
mated the so-called influence parameter by using the low- 
density limit of the direct correlation function (Mayer 
function of the pair potential) to get the density pro- 
file for the drop, and calculated the surface tension from 
the Irving-Kirkwood |94| pressure-tensor expression; a 
monotonic decrease of the surface tension with decreas- 
ing drop radius was predicted (corresponding to 5 > 
throughout) as obtained by Tolman [T5] thirty years ear- 
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lier. The same was found by Hooper and Nordholm [113] 
with a similar generalized van der Waals approach. Guer- 
meur et al. |115) also employed SGT in a similar way to 
Falls et al. |107j . but using a density dependent influ- 
ence parameter, and computed the surface tension from 
an extended Laplace expression: by contrast, these au- 
thors found a non-monotonic dependence of the surface 
tension which increases from below the value of the pla- 
nar interface as the drop radius is increased, becomes 
larger than the planar value and exhibits a maximum 
at about 10 molecular diameters, then decaying slowly 
to the planar limit, corresponding to a small positive 
<5oo. The same overall behaviour as that observed by 
Guermeur et al. (1 15] has now been found in more re- 
cent studies with variants of the SGT approach (e.g., see 
references [HH EM EM E3JJ ) . 

Further controversy has surrounded attempts to in- 
clude higher-order curvature corrections in the expan- 
sion of the surface tension, i.e., to add terms beyond 
the first-order Tolman correction. Strictly speaking, Tol- 



log [ 7 (i?)/7o 



25 /r 2 [l + 6/r+l/3{5/r) 2 ] 
1 + 25/r [l + S/r + l/3{5/r)' 



does involve higher-order contributions, which after ne- 
glecting the terms 0(5/r) in comparison with unity, and 
treating d as a constant, leads to the compact relation 
Q. Tolman himself did not put a firm reliance on his 
expression when considering very small droplets are con- 
sidered. He questioned two assumptions leading to his 
final expression: firstly, that 6 in Eq. ^ is a constant 
for any drop radius; and secondly, the anticipation of a 
bulk liquid behaviour in the centre of the drop. Interest- 
ingly, Tolman suggested that the thermodynamic con- 
cepts should be replaced by "a more detailed molecular 
mechanics" treatment for very small droplets. 

With the goal of generalizing the description to highly 
curved interfaces it is tempting to extend Tolman's the- 
ory to higher order with a formal expansion of the sur- 
face tension in powers of curvature. Helfrich [21j intro- 
duced such an expansion for the surface tension of general 
curved surfaces to second order in the curvature, which 
for a spherical interface can be expressed as 



l(R) = 7oo + 2kC - + (2k + 27?) 
H 



R 2 



(9) 



where Co is the so-called spontaneous curvature, n is the 
rigidity constant of bending, and k is the rigidity constant 
associated with the Gaussian curvature (which is 1/R 2 in 
the case of a sphere) characterizing the energy penalty for 
topological changes of the surface. The original expan- 
sion of Helfrich [3T] is a general form of a second-order 
surface free energy and its derivation was motivated by 
the ultimate goal of describing the elasticity of lipid bi- 
layers that make up cell membranes. It is wildly recog- 
nized as the basic formalism for the description of the 
mechanical behaviour of biomembranes and liquid crys- 
talline phases. Clearly, by including the second-order 
term in the Helfrich expansion one takes the step from 
pure thermodynamics (Tolman's approach) to the the- 
ory of elasticity: while in the Tolman-Gibbs concept 7 
is viewed as an excess (over the respective bulk phases) 
interfacial free energy per unit area, it is the force acting 
against the distortion of the surface in the phenomeno- 
logical Helfrich approach. 

Fisher and Wortis |lllj have used a curvature expan- 
sion of the density and chemical potential with a Landau 
free energy (of square gradient form) to examine the Tol- 
man length. Using an Ising-like model they show that 
<5 = is a general result for models characterized by 
a symmetrical order parameter (density) profile, a con- 
clusion also arrived to by Rowlinson [128] , In the more 
general case of an asymmetric fluid system treated at the 
van der Waals level S ~ — 0.02er, which complements the 
exact results for the penetrable-sphere model (for which 
Hemingway et al. [108] found the exact result 6 = —a/2 
in the zero temperature limit) as the Landau approach is 
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a mean-field theory applicable in the vicinity of the criti- 
cal point. For another typical application of the Helfrich 
curvature expansion within a density functional theory 
the reader is directed to work of Romero-Rochin et al. 

Blokhuis and co-workers |127) 1137) 1149) have exam- 
ined the thermodynamic properties of curved interfaces 
with curvature expansions of the free energy in a series 
of enlightening papers. Making use of the Helfrich cur- 
vature expansions van Giessen et al. |137j found negative 
Tolman lengths for liquid drops at all temperatures, and 
demonstrated that the sign of the Tolman length is very 
sensitive to the details of the free energy in a series of in- 
teresting papers. A negative value of the Tolman length 
has also been predicted on the basis of purely thermo- 
dynamic expansions by Bartell |142j and by Blokhuis 
and Kuipers [149] , the former proposing a simple re- 
lation between the Tolman length and the isothermal 
compressibility ki of the liquid at two-phase coexistence, 
5 « -At/7. 

Controversies associated with the use of curvature ex- 
pansions, and in particular the relevance of the second- 
order correction to the surface tension, follow from the 
fact that the second-order term is proportional to the 
area of the interface and does not therefore contribute 
to the overall free energy, i.e., it just leads to a shift 
of the thermodynamic potential and cannot thus play 
any role in the restoring force acting against the surface 
distortion. There is also evidence of a non-analyticity 
in the free-energy curvature expansion |U 11011 11281 11351 
11441 11451 11551 1156j which suggest that the expansion of 
the surface tension in R^ 1 is generally inappropriate be- 
yond the leading-order term. Studies of fluids in contact 
with hard spherical substrates lead to the conclusion that 
there is a non-analytical contribution of the In R form 
[1351 11441 11451 1156] ; though such a system is clearly not 
the same as a free liquid drop one may expect a cur- 
vature dependence of this type for particles with long- 
ranged interactions, particularly in the vicinity of the 
critical point, but this would be very difficult to identify 
in practice. 

It was recognized early on that the most promising 
route to understanding the intricacies of curved surfaces, 
and liquid drops in particular, would involve a rigorous 
microscopic statistical mechanical treatment. Classical 
density functional theory has amply proved to be a pow- 
erful tool for the description of the intcrfacial properties 
of fluids [25], and is therefore a particularly appropriate 
approach. One of the first applications of DFT for liq- 
uid drops was by Lee et al. |116j who employed a mean- 
field perturbation theory in the canonical ensemble with 
a local density approximation (LDA) for of the hard-core 
reference term (MF-DFT). As will be re-enforced later 
in our paper, the advantage of the canonical ensemble 
is that one can study "stable" equilibrium droplets to 
provide the thermodynamic and structural properties of 
the system [1101 1112j . Lee et al. [116) evaluated the in- 
terfacial tension of the drop using a combination of the 



Laplace and Tolman relations with the pressure tensor 
at the centre of the drop as the corresponding value of 
the internal liquid pressure (obtained locally by identi- 
fying the tangential component of the pressure as the 
negative of the grand potential). This approach leads 
to a monotonically decreasing dependence of the surface 
tension from the planar limit with increasing curvature, 
and correspondingly a positive Tolman length; though 
the extrapolated value for the planar limit of the func- 
tion 5(R) appears to tend to zero, the corresponding er- 
ror bar is large. In the subsequent work of Talanquer 
and Oxtoby [131] with a similar MF-DFT approach, a 
small negative value of the Tolman length was obtained 
by extrapolation, but again a near monotonic decrease of 
the surface tension with curvature was found. Both the 
Lee et al. |116) and Talanquer and Oxtoby [131] stud- 
ies were carried out in the canonical ensemble as it is 
then straightforward to stabilize the drop in a finite-sized 
system. By contrast, Oxtoby and Evans |118) studied 
the nucleation of liquid drops with MF-DFT in an open 
system (grand canonical ensemble), and determined the 
barrier of nucleation from the maximum in the grand 
potential as a function of the supersaturation (drop ra- 
dius). The predictions of the MF-DFT for the barrier in 
the grand potential were compared with those obtained 
with the classical nucleation theory (CNT) (which re- 
quires the planar vapour-liquid tension as input): the 
barrier height obtained from MF-DFT was lower than 
the value obtained from CNT in the case of small drops, 
and was seen to increase above it as the drop size was 
increased. As there is a direct link between the barrier 
in the grand potential (work of drop formation) and the 
surface tension, the finding of Oxtoby and Evans [118] 
implies that the tension of the drop rises above that of 
the planar limit (which would thus be consistent with 
a maximum in the surface tension and a negative Tol- 
man length). In a subsequent paper Zeng and Oxtoby 
|123j have extended the treatment for the more realistic 
Lennard- Jones potential, and good agreement is found 
for the condensation nucleation rates of nonane. 

For the sake of a mathematical convenience, Oxtoby 
and Evans |118j applied the Sullivan hard-core Yukawa 
model [105] . The advantage of using such a model is that 
the Euler-Lagrange equation corresponding to the mini- 
mization of the grand potential functional can be written 
down in the form of a differential equation that is easier to 
solve than the integral equation obtained from the stan- 
dard variational approach. In contrast to the original Sul- 
livan study of a planar interface, however, the boundary 
conditions for the spherical geometry are much less ob- 
vious for the liquid phase. The Sullivan MF-DFT model 
was also adopted by Hadjiagapiou |129j and following a 
mechanical (pressure-tensor) route a non-monotonic de- 
pendence of the surface tension as a function of the drop 
radius was found. This is consistent with the findings 
of Oxtoby and Evans [118] . but in contradiction with 
those of Lee et al. [116] . However, the surface tension 
reported by Hadjiagapiou [129j is higher than that of 



the planar surface over the whole range of radii consid- 
ered, and S(R) — R e — R s < is found to decay almost 
linearly with increasing drop radius which is rather sur- 
prising particularly in view of its magnitude (5 ~ 10 for 
R = 50cr). 

Instead of using a mechanical approach that suffers 
form the ambiguity of the definition of the pressure ten- 
sor, Koga et al. |136) have undertaken a very clear and 
thorough DFT study of both the Lcnnard- Jones and 
Yukawa models within the LDA, analyzing the surface 
properties of the liquid drop on the basis of the Gibb- 
sian thermodynamic theory of capillarity. In qualitative 
agreement with the work of Hadjiagapiou [129] . Koga et 
al. obtained a non-monotonic behaviour for the surface 
tension and a negative Tolman length, but the functional 
dependence found for 5(R) is very different: there is a 
rapid increase in S on decreasing R for R < 10a with 
a change in sign from negative to positive suggesting a 
rapid decrease of the surface tension below its limiting 
planar value for drop diameters corresponding to a few 
molecular diameters. For R > 10cr the length S(R) de- 
cays very slowly to its asymptotic value (Tolman length) 
reaching a magnitude of about one tenth of the molecular 
diameter. However, the authors still consider the sign of 
the Tolman length to be elusive and suggest this merits 
further investigation. 

All of the DFT studies mentioned thus far follow a lo- 
cal treatment of the reference free energy, neglecting the 
short-range correlations in density an subsequent inho- 
mogeneities which may be important in the case of small 
drops. A weighted density approximation (WDA) can be 
used to incorporate these correlations in the free energy 
functional. The first to employ this type of non-local 
DFT were Bykov and Zeng [Mil EH DIE] who used the 
WDA-DFT of Tarazona [161] combined with the gener- 
alized formula for the surface tension and Tolman length 
of Blokhuis and Bedeaux [127] . A non-monotonic curva- 
ture dependence and negative Tolman length were found 
by Bykov and Zeng, though rather surprisingly the dif- 
ference between the WDA and the LDA treatment was 
rather small. In more recent work Li and Wu [154) have 
used a non-local DFT, the fundamental measure theory 
(FMT) of Rosenfeld [162] . to treat the hard-core refer- 
ence perturbation term, together with a quadratic ex- 
pansion of the attractive contribution to the free energy 
where the direct correlation function is described with 
the mean-spherical approximation (MSA). In contrast to 
the findings of a number of the other DFT studies (cf. 
Refs. d2H [US] EH ES]), Li and Wu [HI] reported a 
monotonic decrease in the surface tension with increas- 
ing curvature, a feature that is consistent with the early 
LDA-DFT work of Lee et al. [116J. However, unlike 
Lee et al., the Tolman length calculated by Li and Wu 
is negative, which appears to be inconsistent with the 
behaviour observed for the curvature dependence of the 
surface tension. In more recent calculations with a simi- 
lar FMT-DFT El [73J EE] a non-monotonic dependence 
of the surface tension with curvature was found, and the 



Tolman length was calculated to be small but negative. 

It is useful at this stage to summarize the rather mud- 
dled state of play of the work involving DFT calcula- 
tions: the non-monotonic behaviour and weak maximum 
in the surface tension observed with varying drop radius 
in the latest FMT-DFT studies El El EH] are in line 
with the findings of much of the other work employing 
the extension of the Sullivan model to a spherical ge- 
ometry within LDA-DFT \TjM EH E3E] and non-local 
WDA-DFT P3SD331CEIH] approaches, and with the lat- 
est simulation data |681 1721 173) , but are in contradiction 
with the results of the DFT studies by Lee et al. |116j . Li 
and Wu EM], Zhou et al. QH], and Corti et al. EST]. In 
the case of bubbles, the FMT-DFT calculations of Binder 
and co-workers El EH] lead to the expected monotonic 
decay of the surface tension from the planar limit as the 
radius of the bubble is decreased; this corresponds to a 
negative Tolman length as obtained for drops. We should 
note that Binder and co-workers perform their variational 
analysis in the grand canonical ensemble which involves 
locating a saddle point in the free-energy surface. As 
Oxtoby and Evans |118) have pointed out it requires a 
specific numerical procedure. An analysis in the canoni- 
cal ensemble is simpler as this involves the minimization 
of the free-energy functional [116] . We will discuss full 
details of our analysis of the curvature dependence of the 
interfacial properties of both drops and bubbles with a 
FMT-DFT treatment in the canonical ensemble in later 
sections of this paper. 

It is apparent that the final conclusions of the large 
body of theoretical work on the curvature dependence of 
the surface tension and the sign and magnitude of the 
Tolman length is still a matter of controversy. This is 
also true of the conclusions drawn from the correspond- 
ing simulation studies. In our paper we return to the 
main question of the curvature dependence (monotonic 
or non-monotonic) of the surface tension, the sign of the 
Tolman length, and the applicability of the Tolman equa- 
tion. We will show that it is not just a matter of choos- 
ing the appropriate simulation methodology (pressure- 
tensor route, free-energy calculation etc.) or theoretical 
approach (SGT, LDA-DFT, FMT-DFT etc.), but that 
the specific analysis of the interfacial properties includ- 
ing the density profile and excess free energy is of key 
importance. We start by making some general obser- 
vations regarding purely mechanical approaches (Section 
2), where we show how the surface tension and Tolman 
length can be represented with a classical Newtonian pic- 
ture. In Section 3 we revisit the main developments of 
the Gibbsian theory for the thermodynamics of spheri- 
cal interfaces, and discuss the key features of the Tol- 
man approach. A novel generic expansion of the ther- 
modynamic relations in terms of the curvature of the 
drop is also developed in this section. The more de- 
tailed molecular-level statistical mechanical approaches 
are discussed in Section 4, including the pressure-tensor 
(mechanical) routes to the surface properties and density 
functional theories. Specifics of a non-local (FMT) ap- 
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proach in the canonical ensemble are then described. The 
numerical calculations with our FMT-DFT for drops and 
bubbles are presented in Section 6, and a detailed analy- 
sis of the theoretical results is made from the macroscopic 
mechanical and thermodynamic perspectives in order to 
assess the applicability of the various routes to the inter- 
facial properties of systems with curved surfaces. 



II. MECHANICAL APPROACH 

A mechanical treatment of interfacial properties dates 
back to the beginning of nineteenth century when an un- 
derstanding of the behaviour of matter relied entirely on 
Newtonian classical mechanics. It was therefore natu- 
ral to explain phenomena such as a capillary rise from a 
mechanical perspective based on the assumption of a uni- 
form distribution of molecules interacting via strong and 
short-ranged (compared to the gravity) attractive forces. 
This followed from the observation that the height and 
the curvature of the meniscus of a liquid in a small capil- 
lary is independent of a thickness of the material making 
up its walls. A crude mechanical treatment of matter, 
though unsuitable for a description of interfacial prop- 
erties of very small droplets, can still provide some in- 
sight on the link between intermolecular forces and the 
macroscopic properties of liquids based on an exclusive 
application of Newtonian physics. 

The existence of a surface tension at a liquid interface 
was recognized in the earliest studies of interfacial phe- 
nomena. In the following development of a entirely me- 
chanical expression for the surface tension we generalize 
the formal approach of Laplace [2] , Dupre [164j , Maxwell 
[165] , Rayleigh [91] and others (as exposed so beautifully 
by Rowlinson and Widom, see [3] and references therein) 
for the work associated with the separation of two planar 
liquid surfaces to form a spherical cavity. This allows us 
to obtain purely mechanical expressions for the surface 
tension and the Tolman length of liquid drops. A molec- 
ular concept of the surface tension can be established on 
the basis of a mechanical equilibrium condition assuming 
the existence of pairwise additive attractive interactions 
u(r) between molecules where the integral 



$ = -p 



drw(r) = 2np / drr 2 u(r) , (10) 
Jo 



is taken to express the mean-field cohesive energy per 
particle. Here, we further assume that u(r) is only a func- 
tion of the radial distance and negligible beyond a certain 
cut-off distance d (so that u(r) = for r > d) which is 
small compared to the size of the system. Furthermore 
one assumes that the number density p is constant, i.e., 
that correlations between particles are neglected (mean- 
field approximation). The latter requires, in particular, 
that the integral in Eq. ( 10 ) is taken in the isotropic part 
of the liquid, at least within the range of u{r). A su- 
perficial particle (one at the interface between the liquid 
and its vapour, the density of which is neglected in our 



current development) lacks some portion of the cohesive 
energy compared to a particle in the interior due to a 
lower number of neighbours and is therefore in a state 
of higher potential energy. This, in turn, means that in 
the absence of an external field the liquid will strive to 
minimize its surface area. The surface tension 7 can then 
be defined as a work that has to be done to increase the 
area of a liquid surface by unit area, 



5W = jSA , 



(11) 



or, alternatively, as the restoring force per unit length 
acting against an increase in surface area. 

The radius of a mechanically stable liquid drop can 
directly be determined from the principle of virtual work. 
Let pi and p v be the (scalar) pressures of the (interior) 
liquid and the (exterior) vapour phases, respectively. The 
work necessary to bring about the change in volume due 
to an infinitesimal isotropic expansion is 



5W V = (p v -pi)A6R, 



(12) 



where A = AttR 2 is the area of the unperturbed surface 
and SR represents the displacement of the surface to- 
wards the vapour phase. Such an expansion of the drop 
leads also to an increase in the surface area by an amount 
5 A = 8ttR5R producing a corresponding surface contri- 
bution to the work 8Wa ■ The total work due to the 
virtual volume expansion in the drop radius is therefore 
given by 

SW = SW v + SWa = AttR 2 { Pv ~pi)5R + %hR 1 5R. (13) 

In (mechanical) equilibrium the work has to be zero from 
which one immediately obtains the Laplace relation ([!]). 

By employing a simple static molecular model a link 
between the surface tension and the intermolecular forces 
can be made. To this end, we calculate the work required 
to separate a drop of liquid of radius R from a bulk liquid 
(i.e., the formation of a vacuum cavity of radius R + d 
in a uniform liquid with a liquid drop of radius R in its 
center) by calculating 



SW= / F(t)d£, 



(14) 



where F(£) is the force between two concentric spherical 
surfaces a distance i apart (see Fig 1). In the following 
development we assume d <C R- This type of approach 
has also been used by Fowler [163] to represent the sur- 
face tension of a perfectly sharp vapour-liquid interface. 
Now we set out to calculate the work needed to unbind 
the molecules in outer layer from the drop. The radial 
(and the only nonzero) component of the force between 
the drop and the molecules that are at a distance t from 
the drop surface (see Fig 1) can be obtained from 



F r (£) =-4irp 2 / dr(r + R) 



r+2R 



dsf(s) 



(15) 
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dip, 
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du(s) 



is the force between two molecules 



where f{s) - ds 
a distance s apart, cos#/(s) is the projection of the force 
in the radial direction, and cj) is the azimuthal angle in 
the usual spherical coordinate system. 




FIG. 1. Sketch illustrating the variables for the calculation of 
the surface tension 'y(R) of a drop with a static mechanical 
approach. 



After integration over the angular variables one obtains 
the following expression: 



F r (l) 

Att 2 p 2 



dr(r + Rf 



r+2R 



dsf(s) 



(r 2 + 2rR + s 
4(r + R) 2 



2\21 



f-d rd 

J dr J dsf(s) [r 2 s + 2rsR + 2sR 2 - s^6) 



We proceed by expressing the force as an expansion to 
leading order in the curvature (\/R) about the planar 
limit: 



'"' {Fo ~ r Fi ~°{r:^ ' - (1Tl 



where the reference planar term is 

rd i-d 



dr / dssu(s) 

= 21 [ dssu(s) 

f d d f d 
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dr 
drru{r) 



d 

2 I drr 2 u(r) 
h 



(18) 



and the leading-order curvature correction is 

r-d r-d 

F\ = —2 I drr I dssu(s) 



drr 3 u(r) . 



I \ drru(r) 



(19) 



After substituting Eqs. ( 17 )— ( 19 ) into Eq. (|14j) and 



integrating by parts, one can express the work done in 
creating the cavity around the drop of liquid as 



r-d r-a 

W= F r (£)de = 4TT 2 R 2 p 2 
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d£ 
(20) 



where 



W = 2 d££ drru{r) - 2 1 dl I drr 2 u(r) 
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and 



d r d 

•2 



Wi = / dtt 
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drru(r) 



dl I dr 3 ru(r) 
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dll 4 u(l) . 



(22) 



The expression Wo for the planar limit was already know 
to Laplace [5] and Dupre |164j , but to our knowledge the 
first-order curvature correction W\ has not been devel- 
oped in this manner before. The surface tension corre- 
sponds to the work per unit area, j(R) = , since two 
surfaces with areas AttR 2 and 4tt(R + d) 2 « AirR 2 have 
been created we have 



j(R) 



where 
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is the surface tension of a planar interface, and the coeffi- 
cient proportional to the first-order curvature correction, 



3j^d££ 3 u(£) 



(25) 



is effectively a mechanical representation of the Tolman 
length, which will be defined on thermodynamic grounds 
and discussed in detail in the next section. The expres- 
sion for the planar contribution is of course identical to 
the one obtained when two planar liquid surfaces are sep- 
arated from each other (e.g., see the derivation by Rowl- 
inson and Widom (3J). Our expression for the Tolman 
length is different from that obtained by Kirkwood and 
Buff [53] or by Schofield and Henderson [105]. The Tol- 
man length in our Eq. (20) is expressed as the ratio of 
the fourth and third moments of the pair potential en- 
ergy. Following a pressure-tensor route for the surface 
tension and the surface of tension, Kirkwood and Buff 
[93] obtained a relation for the Tolman length which is 
proportional to the ratio of the fifth and fourth moments 
of the corresponding pair virial (force). Our expression 
(20) yields a Tolman length which is always negative, 
while that obtained by Kirkwood and Buff is of roughly 
the same magnitude but is always positive. For a square- 
well potential of range 1.5a we find <5 ~ — 0.4<r, where a 
is the hard-core diameter, and in the case of the (full) 
Lennard- Jones potential one obtains 6 = — |<7 (where 
the integration is carried out from a up to infinity in 
both cases). We shall return to the issue of the sign of 
the Tolman length later in the discussion. 



III. THERMODYNAMIC APPROACH 

It is evident that the properties of liquids are highly de- 
pendent on temperature. This feature is ignored within a 
purely mechanical perspective, such as the one described 
in the previous section, where the molecules are presumed 
to be at rest in positions of minimum potential energy. 
By the last quarter of the nineteenth century, Boltzmann 
[166) had developed his kinetic theory of gases, a molecu- 
lar theory based on Newtonian mechanics, revealing that 
the temperature of the system is related to the mean- 
square velocity of the chaotic motion of the particles. As 
a result, a dynamic, rather than a static view of molec- 
ular systems began to be accepted. Boltzmann's theory 
provided a dramatically new insight into the behaviour 
of fluids but remained essentially inapplicable until the 
advent of computers and the development of numeri- 
cal molecular dynamics techniques. It is now fully ac- 
cepted that mechanics itself describes only a part of the 
physics of fluids, that directing a system towards its en- 
ergetic minimum, but not its counterpart which demands 
a maximization its entropy. Within a strict thermody- 
namic treatment of the interfacial properties of liquids 
one abandons the molecular picture, describing matter 
as structureless, and providing instead general relations 



between its macroscopic properties. The thermodynamic 
description of finite systems, such as a drop of liquid 
nucleating in a fluid, is far less obvious and cannot be 
properly handled when the dimension becomes truly mi- 
croscopic. However, the concept of the surface tension 
and the related characteristics of the surface can be rig- 
orously formulated within the Gibbsian thermodynamic 
approach without any restrictions, as will be briefly sum- 
marized in the following section. 



A. Theory of Gibbs 

We consider a one-component system containing a liq- 
uid drop surrounded by its vapour. In section [H] we im- 
plicitly assumed that such a system could be character- 
ized by a parameter R, representing the radius of the 
drop. However, such a division of the system into two 
physical subsystems, one corresponding to the liquid and 
one to the vapour phase, is not evident unless the inter- 
face is perfectly sharp, which is never the case. In order 
to avoid this problem, Gibbs [8HTU] introduced a math- 
ematically rigorous theory where one relies on a formal 
definition of a dividing surface separating the system into 
two hypothetical uniform subsystems. Thus, the volume 
V of the entire system is expressed as a sum of volumes 
of the two subsystems 



V = Vi + v v 



(26) 



with the liquid volume expressed as Vi — ^nR 3 . Any sur- 
face variable is now defined as the excess of the variable 
X of the entire system over the sum of the corresponding 
variables of the two bulk subsystems: 



X s — X — Xi — X v 



(27) 



where Xi and X v are the properties of the liquid and 
vapour systems, respectively, at the same thermody- 
namic conditions as the system of interest. The latter 
condition is the key to the thermodynamic treatment of 
interfacial systems. With this division of space in hand, 
one can define the surface adsorption per unit area as 



T(R) 



1 



A(R) 



(N — Ni — N v ) 



N S {R) 
A(R) 



(28) 



where one refers to a given choice of the radius R which 
also defines the dividing surface. 

In terms of thermodynamics, the mechanical model 
adopted in section [IT] corresponds to an adiabatic pro- 
cess with the internal energy playing the role of the ther- 
modynamic potential being minimized by an appropriate 
compromise between the volume and the surface corre- 
sponding to a force balance (cf. Eq. (13l). If the pro- 



cesses are carried out at fixed temperature, which is both 
experimentally relevant and computationally more con- 
venient, the relevant thermodynamic potential is a free 
energy. For a one-component liquid drop, the choice of 
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R defines the liquid and the vapour volumes, Vi and V v , 
and the total differential of the Helmholtz free energy can 
be expressed as 

dF = -piAVi -pidVi+ SdT + pdN + 7 dA + CdR , (29) 

where pi and p v are the scalar pressures of the uniform 
liquid and vapour systems corresponding to a given equi- 
librium chemical potential //, and C is the conjugate vari- 
able to R. In Eq. (29) the use of the general equilibrium 
conditions T = T\ — T v and p = pi = p v has been made 
implicitly. Eq. (29 ) can be integrated over the whole 



spherical surface at fixed R and T [3] 
F = —piVi — p v V v + jA - 



pN 



and from Eq. (27), the surface free energy can be identi- 
fied as 



F s = jA + pN s 



(31) 



so that when referred to a given dividing surface R, the 
surface tension can be expressed as 



i(R) 



A 



(32) 



In other words the surface tension is the surface free 
energy per unit area, providing there is no net adsorp- 
tion r = 0. In a one-component system the choice cor- 
responding to vanishing adsorption is the equimolar or 
Gibbs dividing surface: R = R e . 

In an open or inhomogeneous system it is more conve- 
nient to use the grand potential defined as the Legendre 
transform f2 = F — pN . In terms of the grand potential 
the surface tension is given by 



j(R) 



A ' 



(33) 



regardless of the choice of the dividing surface. 

As the dividing surface is fixed by convention, the 
free energy for fixed N, V, and T or the grand poten- 
tial for fixed p, V, and T cannot depend on this formal 
choice (nor can any other thermodynamic quantity, such 
as pi , p v , p or T). If, for a given drop, we take the formal 
derivative of Eq. ( 30 ) with respect to R and compare it 
with Eq. ( 29 1 , one obtains a generalized Laplace relation, 



Pi - Pv 



g7(g) 
R 



C 
A 



(34) 



[cf. Eq. Q], with an explicit form for the conjugate 
variable: 



C 



dR 



A. 



(35) 



The dividing surface R s for which C = 0, i.e., the one for 
which the macroscopic Laplace relation [cf. Eq. ([T])] is 
satisfied, is commonly referred to as the surface of ten- 
sion at which the tension acts [3]. This requires that the 



formal derivative of the surface tension with the position 
of the surface to be at an extremum: 



d-y(R) 



dR 



= 0. 



(36) 



R=R, 



The generalized Laplace equation can be re-written as 
[i? 2 7(i?)] = R 2 Ap, and on integrating from the sur- 
face of tension R s to an other dividing surface R, one 
obtains [31IT08]: 



7(g) 



R — R s \ R s + 2R 



R 



3-R.s 



(37) 



(30) From Eq. {37} it follows that ^(R) is at a minimum at 

R s ,l(Rs) differs 



the surface of tension, and that, for R 
from 7(i?) by terms of order 1/R 2 . 



B. Theory of Tolman 

Gibbs' theory for the surface tension is based on a def- 
inition of the dividing surface which is taken to sepa- 
rate the two coexisting phases and to which the surface 
tension and the other superficial quantities are referred. 
There are two useful definitions of the dividing surface: 
the equimolar (or Gibbs) dividing surface R e , defined by 
T(R e ) — 0; and the surface of tension R s , defined by Eq. 
( 36 1 . Tolman |13| extended the general thermodynamic 



theory of Gibbs, exploiting the Gibbs-Duhem relation to 
obtain a thermodynamic expression for the curvature de- 
pendence of the surface tension. More specifically, Tol- 
man expressed the adsorption in terms of the difference 
in the two dividing surfaces S = R e —R s . The adsorption 
relative to the surface of tension can be written in terms 
of the appropriate integrals over the number density pro- 
file p(r) as 



1 

AttR 2 

i-O 



(p(r) - pi)r 2 dr + / (p(r) - p v )r 2 dr 



J [p{r + R s )- Pl ] f^ + ^jdr 



[p(r + R s 



p v ][l + —)dr. 



(38) 



Now, if the equimolar dividing surface is a distance 5 
from the surface of tension, R e = R s + S, the adsorption 
r(i? e ) at i? e can be expressed as 

T(R S + 5)= J [p(r + R s ) - Pl ] + ^) dr (39) 
+ J s [p(r + R s )- Pv ](l + ^jdr = 0. 



After combining Eqs. (38) and (40), one obtains [13] 



Pl - Pv 



= 6 



1P_ 
3R 2 



(40) 
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which relates a microscopic property of a drop character- 
ized by the Tolman length <5, to thermodynamic quanti- 
ties that can be obtained directly from Gibbsian thermo- 
dynamics. Indeed, if the Gibbs adsorption equation at 
constant temperature [which follows from Eqs (29) and 
PI], 



d7 = —Tdp , 



(41) 



is combined with the Gibbs-Duhem relations (dpi — pidp 
for i — l,v, at constant T), one obtains 



related to the surface rigidity. The relevance of the l/R 2 
contribution has been questioned |111( 1128] owing to the 
fact that it corresponds only to a constant term in the 
free energy (in three dimensions) and as a consequence 
cannot contribute to the restoring force acting against 
surface deformations. Instead, there is some evidence 
BDH IS51 HS1 H13 H551 H5SI for the existence of non- 
analytic terms, such as ~ lni?, the existence of which is 
still under discussion. 



d7 



Pi 



-d(pi -p v ). 



(42) 



Substituting for the the pressure difference Ap = pi 
from the Laplace equation one can write 



Pv 



C. Curvature expansion 



dry(R) = - 



Pi - Pv 



-d(2j(R)/R) . 



(43) 



This leads to a differential equation for 7(7") in terms of 
the radial integration variable r: 



d7(r) 
7(r) 



' Pi ~ Pv 



r 



r \Pl - Pv 



-dr . 



(44) 



Using Eq. ( 40 1 and integrating the last expression from 



the plane surface (r — 00) to R one finds 

2 r 



In 



7(i?) 

7oo 



' Pi ~ Pv 



1 



r 



r \pi - p v 



-dr . 



(45) 



which with the help Eq. (|40|) at lowest order, i.e., T — 

(46) 



SAp leads to the Tolman equation: 
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7 (i?)=7oo (1-^) +H.O.T. 



In the considerations leading to Eq. ( 46 1 one assumes 



that S is a constant, so that the Tolman length can be 
expressed in the planar limit, such as 



5 = lim (R e — R s ) — z e 

R s — >oo 



(47) 



where the z e and z now define the the corresponding 
perpendicular distances from the interfacial plane. 

The extension of Eq. ( 46 1 beyond the first-order cor- 
rection in curvature is still a matter of controversy. In 
his derivation, Tolman [13_ obtained terms of order l/R 2 
and l/R 3 , cf. Eq. ([8]), but neglected them and expressed 
doubts about their physical relevance taking into account 
the macroscopic origin of his approach. On the other 
hand, in his phenomenological theory, cf. Eq. ([9]), Hel- 
frich [3T] included a l/R 2 'elastic' contribution which he 



We have already showed in Section |TT] that the Tol- 
man length can be represented within a primitive purely 
mechanical standpoint. In the following sections we will 
show how a statistical mechanical treatment can be used 
to provide a reliable and physically consitent estimate of 
5. Prior this, we propose a procedure for the determi- 
nation of the Tolman length from a purely macroscopic 
thermodynamic basis by assuming the analyticity of ther- 
modynamic quantities in the curvature c = R^ 1 of the 
drop. We start by considering a thermodynamic state 
in a metastable region on the vapour side of the phase 
diagram, i.e., a supersaturated vapour with a chemical 
potential which is slightly higher than the saturation 
value. The chemical potential, density and other thermo- 
dynamic functions of such a system can be characterized 
in terms of the radius of the critical nucleus, R = 1/c, 
and we can thus develop a Taylor expansion about the 
planar limit as 



p(c) = p(0)+p'(0)c+-p"(0)c 2 



(48) 



Pi{c) 



Pl (0)+ p i(0)c+-pii(0)c 2 + - 



i = l,v, (49) 



where ' denotes the derivative with respect to the cur- 
vature, and (0) the reference saturation value. As 
the state is supersaturated with p(c) — p(0) > and 
Pi(c) — Pi(0) > 0, it follows that the sum of the first- 



order terms on the right hand sides of Eqs. (48 ) and (49 ) 
must be positive. 

The free-energy density / = F/V of both phases is 
then expanded up to second-order in density making use 
of Eqs. (48) and (49), and the thermodynamic relation 
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df/dp = fi: 



f(Mc)) = f( Pi (0)) 

1 d^l 
+ 2 dpf 



dl 
dpi 



[ Pi (c) - Pi (0)} 



[ Pi (c) - Pi (0)} 2 + 0((pi(c) - Pi (0)f) 



1 d/z 

2 5^ 



^(cl-ftfO^Of^W-fttO)) 3 ) 



/( Pl (0))+/i(0) 



p;(o) C +-pr(o)c 2 



1 9^ 

2 dp~ 



+ C(c 3 ) (50) 



= f(Pi(o)) + cfi(o)p'M 



1 



9/i 



(^(0)) 2 +^'(0)M(0) 



0{c A 



Likewise, using the formal thermodynamic identity re- 
lating the pressure and chemical potential to the free 
energy, both the liquid and vapour pressure can be ex- 
pressed as the corresponding expansions about their sat- 
uration values up to second order: 

Pi(c) = fJ,(c)pi(c) - f(pi(c)) 

p(0)+p'(0)c+^p"(0)c 2 



Pi (0) + ^(0)c+-^(0)c 



-/( A (0))-c M (0K(0) 



,1 



dp 
dpi 



(p'M) 2 + pUoMo) 



+ 0(c 3 ) 



p(o) Pi (o) - /( Pi (o)) 

+c [p(0) P m + p'(0) Pi (0) - A*(0)^(0)] 

2 



(51) 



/(oW(o) + ^"(o) ft (o) + ^(oK(o) 



1 dp 

2 



(^(o)) 2 -Jp'/(o)m(o) 



p(0)+ CA /(0) Pl (0) 



+ G(c 3 ) 

1 ^ 



m'(o)p' j (o) + -m"(o)p 1 (o) o . hi 



{p'MY 



The pressure difference Ap = pi —p v can thus be obtained 
in compact form as 



Ap = cp'Ap + c 2 

1 ,2 

2dp/\ 



+ 0(c 3 ), 



(52) 



where the explicit dependence on the curvature has been 
dropped bearing in mind that all of the terms that are 



retained corr espond to saturation and Ap — pi — p v . Ex- 
pression ( 52 ) can be compared to the combination of the 
Laplace and the Tolman relations, cf. Eqs. ([I]) and (46 1, 
Ap = 27C — 2jSC 2 . On equating the first-order terms 
we obtain 



p'Ap = 2 7 , 



(53) 



implying that p" = —2'y/(Ap) 2 Ap'. One should note 
that Eq. (53) is consistent with the Laplace relation to 
first order [Ee., for 7(c) = 7(0)], as can be seen by com- 
bining Eq. ( 53 1 with the Gibbs-Duhem equation and 



integrating the resulting differential equation from the 
planar limit to some finite curvature. 

An examination of the second-order terms implies that 



27 
Ap 



Ap' 



Ap 2 dpi 



and using Eq. ( 53 ) we find that 



dp 

dpi 



finally arriving at 



27 
Ap 



Ap' 



Ap 



P_ 
P' 



Ap' 



1 dp , 2 
2d Pv Pv 



2l 



-2 7 <S, (54) 



A PP ' t 



Ap 



Ap' = -2 7< 5 , 



where the terms on the right-hand side are seen to cancel 
implying that (5 = 0. 

One therefore reaches the interesting conclusion that 
despite the fact that Tolman's theory is constructed 
purely on thermodynamical grounds, a purely macro- 
scopic thermodynamic treatment yields a trivial solution 
with a vanishing Tolman length: the result implies that 
in the thermodynamic limit, the Gibbs dividing surface 
corresponds to the surface of tension. One can regard 
this apparent paradox as a consequence of Tolman's the- 
ory only providing relations between thermodynamically 
observable quantities, which although well defined at the 
macroscopic scale, their differences are of microscopic di- 
mensions and so beyond the scope of a thermodynami- 
cal treatment. It is interesting to note that Wortis and 
Fisher [lllj also find that 5 — in their analysis of 
symmetric interfaces with a Landau free energy of the 
square-gradient form; the interfaces are symmetrical by 
construction in our purely thermodynamic curvature ex- 
pansion. It is rather ironic, however, that a non-zero 
and physically reasonable representation of the Tolman 
length is obtained with the purely mechanical treatment 
developed in section II, despite the fact that the inter- 
face is assumed be a sharp symmetrical step (see also the 
discussion in Ref. [156J ) . In order to obtain any useful 
information from the Gibbs-Tolman theory one therefore 
has to adopt a molecular (microscopic) approach. 



IV. STATISTICAL MECHANICAL APPROACH 

Statistical mechanical approaches of inhomogeneous 
systems are generally based on determining the response 
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of the system to changes in the external conditions. In 
contrast with the approaches discussed in Sections |TT| and 
|III| statistical mechanics allows for a microscopic treat- 
ment where molecular-level detail can be taken into ac- 
count in a formal manner. 

In this section we outline the common statistical me- 
chanical routes for inhomogeneous systems. These routes 
are not independent and it is important to highlight the 
important interrelationships. One approach (usually re- 
ferred to as the 'mechanical' route) relies on the mechan- 
ical definition of the surface tension as the stress trans- 
mitted across a strip of unit width normal to an inter- 
face. This leads to an expression for the surface tension 
in terms of components of the microscopic stress tensor 
(negative of the pressure tensor) , and can thus be viewed 
as microscopic-level description of the theory of elastic- 
ity. A second approach, the so-called 'virial' route leads 
to an expression for the surface tension which is based on 
the isochoric-isothermal change in free energy due to an 
increase in the interfacial surface by unit area. A third, 
the 'compressibility' route, relies on a calculation of the 
change in free energy arising from an increase in surface 
area caused by density fluctuations. Finally, we present 
a 'thermodynamic' route which allows for the determina- 
tion of the surface tension directly from Gibbsian ther- 
modynamics as presented in Section III A With this dis- 
parate variety of methodologies that are at our disposal 
for a statistical mechanical description of interfacial sys- 
tems it is therefore not altogether surprising that their is 
little convergence in the findings for even the most basic 
of properties. 



A. Mechanical (pressure-tensor) route 

Let us consider a spherical drop in a fixed volume V for 
a system of particles interacting via a pairwise potential 
u(rij) (although this assumption is not restrictive) and 
calculate the instantaneous force on the drop. The force 
is related to the flux of linear momentum density through 
the volume, Fy(t) = J v drJ Q (r, t), where J a (r, t) can be 
expressed as 



forces, 



J a (r,i) = V /3 cr Q/3 (r,t) 



(55) 



when no fields are considered, where the the common 
implicit summation notation of Einstein is used, if not 
otherwise stated. The stress tensor a al3 (r, t) incorporates 
the change in momentum due to particles crossing the 
boundary of V, 



f M) = - 



in. 



S(r - r t ) , 



(56) 



where S(r— r, ; ) is the Dirac delta function, and the config- 
urational part of the stress induced by the intermolecular 



V?af(r, t) = - i E V ?(^') ^ r - r >) - ^ r - r ^ 



(57) 

ct"' 9 itself is not given uniquely, but in general can be 
expressed as |109] 

^(M) = ^E^/ dZ^(r-i), (58) 

i jjii y •' Ci i 

where u'(rij) = &u(rij)/drij, for an arbitrary contour Cij 
joining and ij. 

Defining the pressure tensor as the negative of the time 
average of the stress tensor 



p^(r) = -<<7«"(r,t)>, 



a/3. 



(59) 



and noting that the average of the left-hand side of Eq. 
( 55 ) is zero at equilibrium, one obtains the differential 



conservation law 



VV(r) = 0. 



(60) 



in the absence of external fields. It should also be pointed 
out that the substitution of Eqs. (56) and (57) into Eq. 



(60) and the use of the equipartition theorem leads to 



the first equation of the BBGKY hierarchy [3] . 

From Eq. |57| ) it follows that the components of the 
pressure tensor depend on the choice of the contour join- 
ing the two interacting particles; hence there is an infinite 
number of "acceptable" definitions for the pressure ten- 
sor. This ambiguity can thus be attributed to the prob- 
lem of specifying the portion of the intermolecular forces 
that act across an elementary area [94 . As has been 
discussed at length by Rowlinson [2S], the problem with 
the uniqueness of the definition of the pressure tensor is 
just a particular case of the more general problem with 
the local definition of any many-body thermodynamic 
quantity; the only exception is the chemical potential as 
neatly captured by potential distribution theorem [167J. 
Difficulties associated with the local definition of thermo- 
dynamic quantities do not reveal themselves for uniform 
systems, where any possible ambiguities average to zero, 
but they become relevant for systems with broken sym- 
metry. 

Before discussing the repercussions of the arbitrary na- 
ture of the definition of the pressure tensor, it is worth 
noting that the ambiguity may also be understood from 
a different viewpoint. According to Noether's theorem 
[168] , conservation laws are reflections of the continuous 
symmetry of a given system. In our system, we assume 
translational and rotational symmetry of the intermolec- 
ular potential: the former has been used in the derivation 
of Eq. ([5Tb QUI- Thus, Eq. pi may be viewed as a 



consequence of the symmetry of the Hamiltonian of the 
system defining the conserving current (through Stoke's 
theorem). If there are no further constraints set on p a " , 
then a class of third-rank tensors, the 'superpotentials' 
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q a @ s which are antisymmetric in the last two indices, 
generate an infinite number of pressure tensors differing 
by \7 s q a5 P , all satisfying the condition embodied in Eq. 



( 60 1 . One can take advantage of the non-uniqueness in 
the definition of p a ^ to cast the tensor in symmetric form 
which allows for a definition of the angular momentum. 
The non-uniqueness of the closely related quantity, the 
energy-momentum tensor, is a well known issue in field 
theory and, in particular, in general relativity, where the 
search for the local components of the energy-momentum 
tensor is sometimes referred to as "looking for the right 
answer to the wrong question" [1691 1170] . 

In spherical symmetry, the pressure tensor possesses 
two independent components: 



P(r) = P n (r)e r e r + P t (r)(e e eg + e^e^) 



(61) 



where e r , eg, and are the unit basis vectors, and P n 
and P t are the normal and transverse components, re- 



spectively. Upon substitution of Eq. (61) into the con- 



dition of a mechanical stability, Eq. ( 60 1 , one obtains 

m 



dr 



(r l P n (r)] = r*- 1 [(i - 2)P n (r) + 2P t (r)} , (62) 



for all values of i. In particular, for i = one obtains 



-V : = I dr-{P n (r) - P t (r)} 



(63) 



The integration of Eq. ( 62 1 over the interface gives rise to 



expressions for Ap and, using the Laplace relation, it al- 
lows one to determine the ratio 7 s /i? s , however, it would 
be preferable to determine 7 S and R s independently of 
each other; the latter would provide information on the 
curvature dependence of surface tension through the Tol- 
man relation (47 1. To this end, based on a consideration 



of a force acting on a flat radial strip and the moment 
about the centre of the drop, the following expressions 
can be obtained [5]: 

/■oo 

1 (R S )R S = / [ Plv {r)-P t (r)]rdr (64) 
Jo 



and 



^(R S )R 2 = / Mr) - P t (r)] r 2 dr , (65) 
Jo 

where pi v (r) = piQ(R s — r) +p v Q(r — R s ), with O repre- 
senting the Heaviside step function. In the planar limit 
these expressions simplify to [53"] 



[P n (z) - P t {z)]dz 



and 



z s = — j [P n (z)-P t (z)]zdz, (67) 

Too J— oo 



which are Eqs. (4) and (5) repeated here for convenience. 
The first explicit form of the local pressure tensor for 
a planar liquid-vapour interface was proposed by Irv- 
ing and Kirkwood |94| . who also pointed out its inherent 
non-uniqueness. As an appropriate contour joining the 
two interacting particles they chose a straight line, and 
obtained the normal and tangential components of the 
pressure tensor as 

Pl K (z) =k B Tp{z) -\f dr 12 ^u'(r 12 ) (68) 

A J r 12 



/ dap (2) (ri2,z - az 12 ,z + (1 - a)z 12 ) 
Jo 



and 



P t IK (z) =k B Tp(z) -- ( dr 12 X ' 2 ^V (r 12 ) (69) 

4 J r 12 

x / dap( 2 \ri2, z — az 12 , z + (1 — a) 212) • 
Jo 



From Eq. ( 60 1 it immediately follows that the normal 
component is constant for the planar interface P n {z) = p. 
Evidently, this condition has to be satisfied by any pres- 
sure tensor regardless of the choice of contour. Harasima 
[55] subsequently suggested a different, asymmetric path, 
dividing the vector into parallel and normal compo- 
nents with respect to the interface: 



P t H (z) = k B Tp{z) 



(70) 



-] [ dr 12 ^P^u>(r 12 )pW(r 12 ,z,z + z 12 ) 

4 J r 12 

Furthermore, Harasima |98j showed that the integral 



z n P t (z)dz 



(71) 



is invariant to the choice of pressure tensor for n = 0, 
but not for the higher moments. 



If we return back to Equation ( 62 1 , it follows that for 

i = 2, 



Ap = 



R 2 



drr \pi v - P t (r)] 



(72) 



i.e., the first moment of pi v — P t is invariant to the choice 
of contour in the definition of Pt, but not the higher mo- 
ments [109] . Therefore, with a mechanical route based 
on a microscopic definition of the pressure tensor one is 
unable to determine z s or R s uniquely and cannot there- 
fore provide a consistent way of obtaining the curvature 
dependence of the surface tension. In general, j(R s ) and 
5 determined from Eqs. (64 1 and (65) or Eqs. (66) and 



(66) (67) differ from those obtained from the Gibbs-Tolman 



theory. 

Eq. ( 65 ) was originally derived by Buff [97] , who cal- 
culated the work accompanying a differential increase in 
the area of a spherical segment, keeping the dividing sur- 
face constant. In an open system, the associated work 
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can be identified with the change in the grand potential, 
so that Eq. ( 65 ) may be rewritten as 



An = - P t (r)dr 



(73) 



according to which the transverse component of the pres- 
sure tensor plays the role of the grand potential den- 
sity. This expression has been used frequently (e.g., in 
the LDA-DFT studies of Lee et al. |116j ) since the cal- 
culation of P t allows for a determination of all of the 
thermodynamic and interfacial properties, including the 
surface tension of the drop. However, expression ( 73 ) 



involves the second moment of P t and thus depends on 
the choice of pressure tensor. A thermodynamically con- 
sistence grand potential requires relation ( 73 1 to be in- 



variant with respect to the choice of pressure tensor, in 
essence corresponding to a tautological definition of the 
pressure tensor such that its transverse component corre- 
sponds to the negative of the grand potential functional. 
Thus, the problem can be recast as the need of finding 
the grand potential of a given molecular model. In Sec- 
tion IV. C we show that it can be obtained more directly 
using the compressibility route. 



B. Virial route 

A virial route within a statistical mechanical frame- 
work is based on a generalization of the mechanical for- 
mulae for the work needed to deform a system. This 
route provides a definition of surface tension as the 
isothermal-isochoric change in the free energy (or the 
grand potential in an open system) during a formation 
of a unit area of surface. Note that such a treatment is 
often denoted as a thermodynamic definition of surface 
tension, since it stems from the thermodynamic expres- 



dF 
OA 



(74) 



NVT 



However, here we associate 'thermodynamic' with the 
route based on a determination of the free energy of 
the entire system, such as provided by DFT (cf. Sec- 
tion IV. D). On the other hand one should not confuse 
virial approaches with those stemming from mechanical 
expressions based on the forces acting between the par- 
ticles (cf. Section IV. A). 

The canonical partition function in the limit of zero 
external field can be expressed as 



Z(N,V,T) = 



A 3N 
Nl 



n^dr^ exp 



U({Ti}) 



k B T 



(75) 



If the system is perturbed by a transformation r' — 
r + £(r), the partition function of the deformed system 



acquires the form: 



Z'(N,V,T) 



A 3N 
Nl 



dr'l 



drf 



exp 



k H T 



(76) 

The associated change of free energy to first order in £ is 

IN n \ 

(AF) NVT = /-fc^T^V • fa) + J2t( r i)Vi u ) ■ 

(77) 

A link between the mechanical and virial route can be 
made at first order, if we introduce the configurational 
pressure tensor (r) = — (a"P(r,t)) from Equation (57) 
and substitute it to ((77| . After some algebra one obtains 



AF = — I drp a0 (r)V p u a {r), 



(78) 



corresponding to the well known "stress-strain" expres- 
sion from the theory of elasticity. Considering now a class 
of deformations with zero divergence, so that the change 
of the free energy at first order is associated solely with 
the change of are a |109] . the corresponding expression 
([78} leads to Eq. §63§~ 

MacLellan [95] used the virial route to derive the fol- 
lowing expression for the surface tension of the planar 
vapour-liquid interface: 
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dzi / dr i2 



2 d 2 Ul2 



2 d 2 Ui2 



dz 



'12 



dx\ 2 



P (2) (ri,r 2 ), 



(79) 

first derived by Kirkwood and Buff from the me- 
chanical route, cf. Eq. (66). Later, Lekner and Hen- 



derson |104j reduced Eq. (79) to a simpler three-fold 
integral. Notwithstanding the growing complexity of the 
corresponding algebra, one can go beyond first order in 
£ (cf. Ref. [Hj where the so-called test-area method for 
the planar interface was developed; within this treatment 
one can in principle determine terms of arbitrary order 
in the interfacial free energy). Recently, a free-energy 
expansion due to a perturbativc deformation has been 
applied to spherical liquid drops including the higher- 
order terms |72) , where the change in the free energy can 
be expressed as 



AF = {AU) 
1 



<D(AU A ) 



2k B T 
[(AU 3 



[(AU 2 ) - (AU) 2 ] 

3(AU 2 )(AU + 2(AU 3 ) 



(80) 



The first term in the average of the deformation energy 
has been shown [731 I104| to be equivalent to that ob- 
tained from the mechanical route [cf. the Kirkwood-Buff 
expression [93]]. The numerical results of molecular dy- 
namics simulation of Lennard- Jones fluids revealed that 
the second-order terms of the free-energy expansion do 
not contribute in any appreciable way in the case of the 
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planar vapour-liquid interface. This supports the con- 
sistency between the the mechanical and virial routes 
to surface tension for a planar geometry. It turns out, 
however, that the second-order term in Eq. (80 1 be- 



comes comparable in magnitude (but of opposite sign) 
to the leading-order term for small drops 72 . This is a 
clear consequence of the enhanced effect of fluctuations 
in nanoscale drops when compared to the planar vapour- 
liquid interface. Attempts to describe the interfacial be- 
havior of a microscopic drop by means of the mechanical 
route or first-order virial expressions are therefore clearly 
invalidated. A thorough analysis of the specific role of 
energetic fluctuations on the thermodynamic properties 
of small drops and bubbles will be the subject of future 
work. 



C. Compressibility route 

The virial route to surface tension leads to a statistical 
mechanical expression involving the gradient of the inter- 
molecular potential and the pair correlation function. As 
shown in the previous section, the standard stress-strain 
formulae follow from the first-order change in free energy 
due to a deformation of the area and leads to the me- 
chanical expression of Kirkwood and Buff 93 . Triezen- 
burg and Zwanzig |102| obtained an alternative result in 
terms of the one-body density and the direct correlation 
function. This expression can be derived formally as a 
functional Taylor expansion in the intrinsic free energy 
up to second order in the density distortion due to an 
external field [TU2] : 

ST 



S 2 T 



A ^/ d ^Wr) 

+ i/dr/dr'^(r)^ (r ') Mr) ^ (r0 
= J dr5p(r)[p-ip(r)\ 



+ ■ 



kT 
~~2 



dr J dr'8p{r)5p{r') 



S(r - r') 
p(r) 



-c(r,r') 



(81) 



where T represents the intrinsic free-energy functional. 
From Eq. (81 )the 'compressibility' form of the surface 



tension for a planar vapour-liquid interface can be ob- 
tained as [102] 

7 = ifcT J dz 1 p'(z 1 ) J dr 2 p'(z 2 )(r 2 2 -z 2 2 )c(r 12 ,z 1 ,z 2 ) , 

(82) 

where p'(z) = dp/dz denotes the gradient of the density 
profile (which characterizes the compressibility of the sys- 
tem). Schofield [106] has shown that this expression is 
equivalent to the one of Kirkwood and Buff [S3] , cf. Eq. 



(66). 

Hemingway et al. |108j extended the result to spherical 



interfaces where one finds 



7 



dri / dr 2 p'(r 1 )p / (r 2 ) 
o Jo 



(83) 



< / dr 12 r 12 [rl 2 - (r 1 - r 2 ) 2 ]c(ri 2 ,ri,r 2 ) . 

'|ri— r 2 | 

Expression ( [84] ) has been assessed for the penetrable- 
sphere model {3. 108.. where an approximation for 
c(ri,r2) is available at the mean-field level, and a con- 
sistency between the compressibility route and the ther- 
modynamic expressions for 7 and 8 [cf. Eqs. (!]), ( |46[ ), 
and (47 was found. In the zero-temperature limit the 



model is solvable exactly and in this case the value of the 
Tolman length is 8 = —a/2. 

A connection with the mechanical route can again be 
made by expressing Ap(r) in terms of the strain field 
and by introducing the pressure tensor. Schofield and 
Henderson [5] showed that Eq. (81) reduces to 



drp a/3 

- J drAp Q/3 V^ a (r) +0(£ 3 ). 



(84) 



From this expression one can see that the compressibility 
route captures terms up to second order in £ which can 
be interpreted as capillary waves fluctuations, whereas 
in the first-order virial expression, cf. Eq. (78), these 
fluctuations are absent. 



D. Thermodynamic route — DFT 

As we have seen in the previous section, an expan- 
sion of the intrinsic free-energy functional up to second 
order gives rise to an expression for the interfacial ten- 
sion in terms of the direct correlation function. Unfortu- 
nately, good approximations for c( 2 )(ri,r 2 ) are generally 
not forthcoming, and this puts limits on the applicabil- 
ity of the method. On the other hand, accurate and 
well tested approximations for the full free-energy func- 
tional are now available, that enable one to determine 
the thermodynamic properties of the entire (inhomoge- 
neous) system. Following Gibbsian thermodynamics as 
described in Secti on |III A] the surface tension can be ob- 
tained from Eq. ( 33 ) but now expressed specifically in 



terms of the surface of tension R s : 



■PlVl +PvVy 
47T.R? 



(85) 



Here, pi and p v are the scalar pressures of two hypo- 
thetical bulk phases corresponding to a bulk liquid and a 
metastable supersaturated vapour, characterized by the 
temperature T and chemical potential p; the departure 
of the chemical potential from its saturation value is de- 
noted by dp = p — p s > 0. Expression (85) only provides 
a formal relation between "f(R s ) and R s , and in order to 
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obtain a complete solution of the curvature dependence 
of the surface tension an independent route to either is 
required; a knowledge R s would enable the determina- 
tion of the Tolman length. Separate expressions follow 
directly from Eq. ( 85 ) when expressed in terms of a gen- 
eral dividing surface radius R, such that d ^ J ^~ > =0, 
which implies 



R=R S 



f 3Aft \ 1/3 
\2nAp ) 



and 



/3Aft(A P ) 2 



16tt 



1/3 



(87) 



where again Ap = pi — p v , and Afl = £1 + p v V is the 
work associated with the creation of the liquid drop. 

A DFT approach is based on the construction of a func- 
tional of the one-body density which exhibits a minimum 
at equilibrium that can be associated with the thermody- 
namic grand potential |25j . In the absence of an external 
field, the grand potential functional is of the form 

n[p(r);p] = F\p{r)]-p J drp(v) . 

The intrinsic free-energy functional J-"[/>(r)] can be writ- 
ten as a sum of the ideal and excess J-" ex contributions 



T[p(r)} = k b T / drp(r) [logA^r) - l] + F ex [p(r)}, 



(89) 

where A is the de Broglie wavelength. Variations of J-" ex 
with respect to density distribution provide correlation 
functions of arbitrary order. 

For a thermodynamically stable state the second vari- 



ation 
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S£lr)Sp(r') mus t be positive, so that the solution of 
Eq. (88) is stable with respect to small perturbations. 



For a macroscopic two-phase system this scenario is real- 
ized only for a planar (Sp = 0) interface; here, however, 
we are concerned with a drop in a thermodynamically 
metastable state (dp > 0), which is unstable with re- 



spect to a uniform liquid and thus 



< 0. As 



Sp(r)8p(r') 

a consequence a drop placed in an open system with a 
radius which is smaller than the so-called critical radius 
will evaporate, while larger drops will grow in an un- 
bounded manner resulting in the complete condensation 
of the system. In order to stabilize the drop, we con- 
sider a closed system, characterized by a finite number 
N of particles, with the free energy as the thermody- 
namic potential. As an alternative one can consider a 
weak spherically symmetric external field </>ext(?")> even- 
tually taking the limit <p ex t{r) — > 0. The constraint of 
fixing the total number of particles prevents the unlim- 
ited growth of the drop as this would lead to a depletion 
in the vapour phase and thus to a decrease in the un- 
dersaturation. In our current work, the attractive part 



of the intrinsic free-energy functional is approximated as 
a perturbation from a hard-sphere reference fluid at the 
mean-field level: 

F[p(r)} = k B T J drp(r) [logA 3 p(r) - l] + F hs [p(*)} 

+\J**J dr'p(r)p(r') Ua tt(|r - r'|) , (90) 
Our model thus consists of a hard-sphere repulsive core 



(86) giving rise to JF^s in Eq. ( 90 ) and an attractive term rep- 



resented with a truncated Lcnnard- Jones potential ac- 
cording to Weeks-Chandler-Andersen |171j perturbation 
theory: 




^min ^ T <~l V c 

r > r c , 



(91) 



with r m j n = 2 1 / 6 ct, and r c — 2.5a will be taken as the 
cutoff distance in our calculations. 

For the hard-sphere contribution to the free-energy 
functional we employ the approximation proposed by 
Rosenfeld |162j in his fundamental measure theory: 



Jhs[p(r)] = k B T J $({n Q }), 



(92) 



where the free-energy density $({n a }) can be expressed 
in terms of weighted densities defined as 

n a (r)= dr'p(r)'w a (r-r') a = {0,l,2,3,«l,u2} 



(93) 

where ^(r) = 9(f — r), lUaC 1 *) = 6(f — r), w\(r) — 
w 2 (r)/47rf, w (r) = w 2 (r) / '4nf 2 , w v2 (r) = f,S(r - r), 
and w v i(r) = w 1 ,2(r)/47rf; the hard-sphere radius is set 
to r = a/2. An FMT treatment is deemed necessary for 
small drops because of the large density oscillations which 
will be present particularly in the central liquid region at 
the lower temperatures, corresponding to higher density 
states. 



The general expressions given by Eq. ( 93 ) can be sim 



plified significantly for some particular geometries. In 
the case of a perfectly spherical drop, the density varies 
only in the radial dimension r, so that the calculation of 
the integrals in Eq. (93) reduces to a problem of one- 



dimensional quadrature: 



n 3 (r) 
n 2 (r) 



drV [ 

\r—r\ 

o ~ pr-\-r 
27TV 1 

r 



(r - r') 2 ] p(r') 



dr'r' p(r ) 



(94) 
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>(r') 



< 2 (r)=nl 2 (r)=0 



The equilibrium density profile is found by minimizing 
the free-energy functional, Eq. (89), which leads to the 
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following Euler-Lagrange equation: 

k B T log A 3 p(r) + SFh / [ ^ )] + / drV^att(|r-r'|)=0, 
WJ J 



subject to the constraint 

J drp(r) = N . 



accurate equation of state for hard-sphere fluid mixtures 
proposed by Boublfk |173j is used (see Ref. |174j for de- 
tails) . 

All the quantities are expressed in reduced units: r* — 
(95) r/cr, R* = R/a, p* = pa 3 , T* = k B T/e, and 7* = 7 cr 2 /e, 
where a and e are the size and energy parameters of the 
Lennard- Jones potential. 



(96) 



Eqs. (95) and (96) can be solved self-consistently from 
TVexp [cW(r)] 



p(r) 



4n J^ d drr 2 exp [cW (r)] 



(97) 



where is cW(r) = — SFcx ^g^j kBT is the single-particle 
direct correlation function. The latter can be separated 
into repulsive and attractive contributions as 



c (1) w = E c « 1) ( r )+ c ittW, 



(98) 



where c« (r) — -g^- ® w a (r) has the same form as n a (r) 
with p(r) replaced by and 



c attO" > r c) 



ci«(r < r c ) 



2tt 
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2n 
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dr'r'p(r')£(r,r'), (99) 



dr p(r )r £(r, r ) 







r c +r 



dr'r'p(r')t(r,r'), (100) 



with Z(r,r') = fl';_ r , l dr"r"u att (r")/k B T and £(r,r') = 

dr"r"u M (r"). 
In the following section we present the findings of cal- 
culations with our non-local FMT-DFT for liquid drops 
of varying size. A detailed analysis of the density profiles, 
the coexistence densities of the vapour and liquid regions, 
the curvature dependence of the vapour-liquid tension, 
and the Tolman length is undertaken, making appropri- 
ate comparisons with existing work wherever possible. 



V. NUMERICAL RESULTS 

In this section we present an analysis of the interfacial 
properties of small liquid drops surrounded by vapour 
for a one component system. If not stated otherwise, the 
description is obtained with the FMT-DFT described in 
Section IV. D in the canonical ensemble for a Lennard- 
Jones 12-6 WCA potential truncated at r c = 2.5er [cf. 
equation (91 )] by solving Equation (97) using a standard 
Picard iteration method. In order to represent the non- 
local functional for the hard-sphere reference potential, 
a modified version Rosenfeld's FMT based on the highly 



A. Structure of a microscopic drop 

In Figure 2 we present density profiles obtained from 
our non-local mean-field DFT for three temperatures 
T* = 0.7, T* = 1, and T* = 1.2 (corresponding to re- 
duced temperatures: T r = T* /T* = 0.526,0.752, and 
0.902, respectively). The latter temperature is already 
rather close to the critical point of the bulk vapour-liquid 
coexistence, T* — 1.33. One can demonstrate several 
general characteristic features of the structure of a liquid 
drop from these profiles. At a temperature close to the 
triple point, which in this case occurs at T t * ~ 0.6, the 
structure inside the drops exhibit strong undamped os- 
cillations which extend from the surface to the centre of 
the dense liquid. This type of highly correlated structure 
in the dense interior of the drop can not be accurately 
described with a traditional square gradient treatment 
(e.g., the work of Falls et al. [107j . Guermeur et al. |115j . 
and the more recent papers (1301 11321 11381 fT^Tj ) or local 
density functional theories (e.g., the studies by Lee et al. 
[ITS] , Oxtoby and co-workers [TTS1 QH CCS] , Hadjiaga- 
piou [129] . and Koga et al. |136j . Note that it is the 
non-local character of our density functional, Eq. ( 92 ) , 



which enables one to capture this type of fine structure. 
These oscillations can also be observed to a lesser de- 
gree on the liquid side of a planar vapour-liquid interface 
|176j . but in the case of a spherical interface both the 



amplitude and range of the oscillation is significantly en- 
hanced. It is clear that there is a strong inhomogeneity in 
the density of small drops at low temperatures, and any 
assumption of a uniform liquid region would evidently be 
unrealistic. In particular, one should note that the den- 
sity at the centre of the drop depends on the amplitude 
and wave-number of the almost periodic density profile 
at r = 0, and on the size of the drop. At the intermedi- 
ate temperature of T* — 1, the oscillations almost vanish 
and the density profiles become monotonically decaying 
functions. This means that in the temperature interval 
Tp W £ (0.7; 1) there is a crossover between an oscillatory 
and monotonically decaying density profile corresponding 
to a Fisher- Widom line [17511176] : the construction of the 
Fisher- Widom (FW) diagram is beyond the scope of the 
present paper. Another salient feature is that the density 
p(0) in the centre of the drop remains somewhat higher 
than that of the saturated liquid density and is nearly in- 
dependent of the size of the drop up to some value of the 
drop radius where a sudden decrease in p(0) occurs. For 
the highest value of the temperature considered here, we 
observe a very diffuse interface between the drop and the 
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vapour; as a consequence any approximation based on the 
assumption of a sharp interface would presumably lead 
to a quantitatively unreliable results as one approaches 
the critical point. Nonetheless, we stress that a Gibbsian 
thermodynamic treatment involving the mapping of the 
system into two uniform regions separated by a well de- 
fined dividing surface is free of any ambiguity and thus 
fully applicable regardless of the drop size. The density 
in the centre of the drop is clearly rather sensitive to the 
size of the drop; for a drop of intermediate size corre- 
sponding to the system with N — 800 at T* = 1.2 the 
density at the centre of the drop is essentially the same as 
that of bulk system so that the vapour-liquid coexistence 
crosses the bulk binodal curve at this point. 

In order to shed further insight into the structure of 
the liquid drop, we depict in Figure 3 the dependence 
of the liquid density on the drop radius, taken as that 
corresponding to the Gibbs dividing surface R e . Two 
definitions are commonly considered to define the liquid 
density of a microscopic drop. In one, the density at 
the centre of the drop p(0) is frequently interpreted as 
the liquid density in computer simulation studies. In the 
other common choice a thermodynamic definition of the 
liquid density pi(T, p) is taken, i.e., one corresponding to 
that of a hypothetical bulk phases, cf. Eq. (30), with 



the same chemical potential and temperature as the sys- 
tem containing a drop. For relatively large radii, the two 
definitions of the liquid density practically coincide and 
exhibit a monotonic curvature dependence in line with 
that predicted from the Laplace equation (see Figure 3) . 
However, below a drop radius of R* ~ 10cr a striking dif- 
ference between the curvature dependence exhibited by 
these two densities becomes apparent. While pi(T,p) re- 
mains monotonic, p(0) exhibits a maximum and its value 
eventually drops below the bulk saturation density. The 
presence of the maximum reflects a non-monotonic curva- 
ture dependence of the surface tension, as will be shown 
in the subsequent discussion. Two opposing effects thus 
determine p(0) at small R: the linear increase of the cap- 
illary pressure with curvature due to the factor 1/R in 
the Laplace equation; and the decrease of the surface 
tension for small R (the latter is a surface contribution 
~ R 2 that becomes dominant for sufficiently small R) . e 
On the other hand, the thermodynamic definition of the 
liquid density is merely controlled by the value of the 
chemical potential, i.e., by the measure of the extent of 
supersaturation. Following an isotherm from the binodal 
to the spinodal (the limit of thermodynamic stability), 
the critical radius of the metastable drop decreases while 
the corresponding liquid density must increase. 

The difference between the thermodynamic definition 
pi and p(0) is also apparent from an inspection of Figure 
4, where we compare the vapour-liquid phase coexistence 
behaviour for systems of different size, with stabilized 
drop radii ranging from 3<r to 26a. The vapour branches 
of the finite systems are shifted towards higher densities 
in all cases in a manner commensurate to the correspond- 
ing curvature of the drop (cf. the upper panel of Figure 
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FIG. 2. Density profiles of liquid drops of Lennard-Jones fluid 
as predicted with our nonlocal canonical mean-filed FMT- 
DFT calculations (cf. Section IV. D) for three temperatures 
and system sizes. The system size is controlled by fixing the 
number of particles inside a spherical container of fixed radius 
D — 10a. The horizontal lines denote saturation densities 
of the bulk vapour and liquid phases for the corresponding 
temperature. 



22 




R s * 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



FIG. 3. The density in the centre of the liquid Lennard- 
Jones drop p(0) and the liquid density of the corresponding 
hypothetical bulk phase pi (T, p) as a function of the drop 
radius obtained with our canonical mean-field FMT-DFT (cf. 
Section IV. D). The LJ system is at a temperature of T* = 1. 

4). Evidently, the vapour pressure and thus the density of 
a drop must be larger than the saturation pressure, and 
the difference is described by the Kelvin relation If 
the density of the drop is defined thermodynamically as 
in the lower panel of the Figure 4, there is a similar shift 
of the whole of the coexistence envelope to the right. In 
all cases the vapour is supersaturated so that Sp > and 
as a consequence the density of the liquid phase must also 
be higher than the saturation value. When the density of 
a liquid drop is associated with the value p(0) of the den- 
sity profile at the centre the scenario is quite different. 
For a sufficiently small drop, p(0) may decrease below the 
saturated liquid density, as has already been observed in 
Figure 2. For the smallest system (N — 300) shown 
in upper panel Figure 4, the liquid branch crosses the 
binodal at T* ~ 0.95 since at these conditions the drop 
is sufficiently small. We should note that the "critical 
point" of the drop (if one is able to define the instability 
of the drop in this way) is always lower than that of the 
bulk fluid. 
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FIG. 4. The vapour-liquid coexistence phase diagrams in 
the density-temperature plane obtained for a Lennard- Jones 
drops with our canonical mean-field FMT-DFT (cf. Section 
IV. D). In the upper panel, the liquid density is defined as the 
density p(0) in the center of the drop, while in the lower panel, 
the liquid density p(T, p) is determined according to the the- 
ory Gibbs for a corresponding hypothetical bulk phase. The 
phase behaviour of three finite systems is complemented with 
that for the essentially infinite planar vapour-liquid interface. 



B. Surface tension and Tolman's length 



1. Surface tension 

One of the most advantageous features of DFT is that 
once the solution of the Euler-Lagrange equation, Eq. 
(971, for the density profile which minimizes the grand 



potential Q (or free energy F) is known, the surface ten- 
sion of the drop can be obtained directly from Eq. ( 33 ) 



(or (32)), since Q (or F) is a direct output of the theory. 



In this way, a direct thermodynamic route to the de- 
termination of the curvature dependence of the surface 
tension and Tolman length can be followed, without the 
necessity to determine ill-defined local thermodynamic 
functions. We should note however that the use of local 
thermodynamic routes within a DFT treatment has been 
commonplace (e.g., see Refs. |116l 1154] ). Alternatively, a 
knowledge of the equilibrium density profile allows one to 
calculate the Tolman length by making use of the ratio 
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7(-R)/7oo according to the Tolman relation, Eq. (46) 



1.0 



FIG. 5. The deviation of surface tension 7(-R) from the planar 
limit 7oo as a function of the Lennard- Jones drop radius cor- 
responding to the Gibbs dividing surface. The prediction of 
our canonical mean-field FMT-DFT following the thermody- 
namic route (continuous curve) are compared with the simu- 
lation results: test-area deformations in the canonical ensem- 
ble, Sampayo et al. \72 (circles); grand canonical ensemble, 
Schrader et al. |68| (dashed); molecular dynamics simulation 
following the mechanical route, Vrabec et. al. [63] (squares). 
The LJ system is at a temperature of T* = 0.8. 

In Figure 5 we compare our canonical mean-filed FMT- 
DFT (cf. Section IV. D) results for the curvature depen- 
dence of surface tension with the recent simulation data 
obtained from canonical [72] and grand canonical [68] 
simulation. The most important observation that can 
be gleaned from Figure 5 is that according to both our 
DFT and the simulation data the surface tension of the 
LJ drop is characterized by a maximum between R ~ 5<r 
and R ~ 10cr. This finding is clearly inconsistent with 
studies reporting a monotonic curvature dependence of 
the surface tension. In particular, one should single out 
any mechanical treatment based on a computation of the 
pressure-tensor components (e.g., Eq. (64)) including the 
majority of the simulation studies which follow the me- 
chanical approach originally presented by Thompson et 
al. [5U]. As can be seen in Figure 5 surface tension ob- 
tained by Vrabec et al. [63 from a very thorough molec- 
ular dynamics study following the mechanical route is 
in contradiction with the latest simulation data (and our 
DFT prediction) : not only is the curvature dependence of 
the tension seen to be monotonic throughout, but the nu- 
merical values are up to 25% lower than the more recent 
calculations following a thermodynamic route 68, 72, 73J. 
As was mentioned in the introduction, theoretical square- 
gradient and density functional theories have suggested 



both monotonic (e.g., |107U116lll54j ) and non-monotonic 
(e.g., |115l ll36, 149:) curvature dependencies for the sur- 
face tension. It is however surprising to note that the 
most sophisticated study to date [T53] (FMT-DFT that 
goes beyond the mean-field approximation for the attrac- 
tive contribution) suggests a monotonic dependence of 
the surface tension with curvature; this is inconsistent 
with our DFT results and the latest simulation studies 
(most likely, this is due to the use of classical nucleation 
theory as a connection to the surface tension which is 
known to break down for small drops [118] ). The fact 
that j(R) > 7oo over the whole region of R beyond the 
maximum suggests a negative Tolman length, an observa- 
tion which is again in conflict with the predictions from a 
mechanical treatment. On the other hand, for drop radii 
below the maximum in the surface tension, the surface 
tension steeply decreases below its planar limit. It can 
thus be supposed that the radius corresponding to the 
maximum of j(R) sets a limit to the validity of the Tol- 
man relation which in view of the definition of Soo , leads 
to a monotonic behaviour for the curvature dependence 
of j{R). This can be assessed with a direct calculation 
of the Tolman length using Eq. (47) with the equimo- 



lar surface (obtained directly from the density profile) 
and the surface of tension (cf. Eq. (86)). A compari- 
son of the curvature dependence of the surface tension 
as obtained directly from the Gibbs- Tolman theory (cf. 
Eq. (46)) is made in Figure 6 for the Lennard- Jones 
drops at a temperature of T* = 1. It is seen that the 
difference between the two approaches is nearly indis- 
tinguishable beyond R > 10a, i.e., almost over the whole 
range of radii where j(R) is monotonic. For smaller drops 
with radii below the maximum, the descriptions with the 
full DFT and Gibbs- Tolman approach start to deviate 
dramatically, with an increase in the respective absolute 
slopes but in opposite directions. On the one hand, this 
supports the consistency of a thermodynamic treatment 
for large drops; on the other hand, such an analysis high- 
lights the limit of validity of the Gibbs- Tolman treatment 
which is a macroscopic thermodynamic approach. It is 
important to reiterate that our DFT predictions are in 
a qualitative disagreement with the results based on a 
mechanical approach (i.e., one relying on the pressure 
tensor) , where a monotonically decreasing dependence of 
the surface tension with drop size is obtained (the latter 
corresponding to a positive Tolman length). 

A consequence of the non-monotonic behaviour ob- 
tained for the surface tension is that the assumptions 
leading to the derivation of the Tolman equation must fail 
when the radius of the drop is of order of the range of the 
intermolecular potential. One possible route beyond the 
Tolman equation is to extend the curvature correction 
to the planar surface tension by including higher order 
terms. Alternatively, one can relax the assumption of a 
constant value of S in Tolman equation [1361 11381 1150j , 
and determine the curvature dependent S(R) from Eq. 
(40). The curvature dependence of 5(R) obtained from 



our FMT-DFT (Eqs. (46) and (86)) is displayed in Fig- 
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representing the curvature dependent surface tension for 
use in extended nucleation theories. 



100 120 



FIG. 6. The surface tension of a Lennard-Jones drop as a 
function of drop radius corresponding to the surface of ten- 
sion. The continuous curve denotes the results from a direct 
determination of the surface tension with our canonical mean- 
field FMT-DFT (see Section IV.D), while the dotted curve is 
obtained from the Tolman equation (cf. Eq. (46 1). The 



planar value of the surface tension for planar vapour-liquid 
interface is indicated by the dashed line. The LJ system is at 
a temperature of T* — 1.0. 



ure 7 for truncated Lennard-Jones drops of various size 
at T* — 0.8, for which the Tolman length is determined 
to be S = lim^-j.oo S(R) = — 0.0708er; this value is com- 
pletely in line with the FMT-DFT estimates of Block et 
al. [72]. We observe a steep increase of S(R) at small 
values of R, and an analysis of the the data suggests 
a dependence S(R) = + a/R 2 , which indicates that 
S(R) ~ ^oo for R > lOcr. The lack of a term linear in 1/R 
in S(R) supports the view of Rowlinson |128j that terms 
in 1/R 2 should not contribute to the surface tension of a 
fluid as they do not give rise to a restoring force on de- 
forming the interface; the term in 1/R 3 (corresponding 
to terms in 1/R 2 for S) would of course contribute to the 
surface tension. One should note that terms in h\R/R 2 
for the surface tension have been attributed to the long- 
ranged potentials in the studies of wetting on spherical 
substrates |135l 11451 1155] ; such a logarithmic dependence 
in ~f(R) has not been identified from our FMT-DFT cal- 
culations for the free drops of particles interacting via the 
truncated LJ potential. When our simple quadratic cur- 
vature dependence for delta(R) is introduced in Eq. ^ 
and integrated, the resulting surface tension is in remark- 
ably good agreement with the values obtained from our 
FMT-DFT with the direct thermodynamic route over the 
whole range of radii. This numerical analysis should not, 
however, be taken as an extension of the original theory 
of Tolman, as one can not establish the physical relevance 
of the correction term and, in particular, one is unable 
to predict the value of the constant a. Nevertheless, em- 
pirical approaches of this type could be useful in, e.g., 
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FIG. 7. Upper panel: The curvature dependence of S(R) = 
i? e — R 3 of Lennard-Jones drops obtained from our canonical 
mean-field FMT-DFT with Eq. pOl. Lower panel: Surface 



tension as a function of radius R e corresponding to the Gibbs 
dividing surface. The symbols represent calculations with the 
Gibbs theory [cf. Eq. ([85], the dashed line to the Tolman 



equation ( 46 1 and the full curve to the modified Tolman equa- 
tion with VI [R*) = 5* + a/R* 2 and a = 1.52733. The value 
of the planar limit of the surface tension is denoted by the 
dotted line. The LJ system is at a temperature of T* = 0.8. 

Throughout our computations we have considered a 
molecular model with Lennard-Jones attractive interac- 
tions truncated at a distance r* — 2.5 from the centre of 
the particle. One may ask how the range of the attrac- 
tive forces affects the interfacial properties of small drops 
of liquid. In Figure 8 we plot the curvature dependence 
of surface tension of drops of LJ particles with differ- 
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ent cutoff distances that are frequently used in simula- 
tion studies. The corresponding planar values of the bulk 
vapour- liquid surface tension increase with r c : a larger 
cutoff implies stronger cohesion and thus a higher value 
of the surface tension, as is also apparent from Figure 8. 
Qualitatively, however, the non-monotonic behaviour of 
the surface tension with curvature remains unchanged. 
It is perhaps just worth noting, however, that the ra- 
dius below which 7(i?) < 700 (indicated by the arrows 
in Figure 8) increases with increasing r c . We recall that 
such a crossover occurs when the surface effects begin to 
dominate the forces in the interior volume, i.e., when no 
"bulk" fluid region can be assigned inside the drop. In 
this case even particles in the centre of the drop "feel" 
the interface, a scenario that becomes increasingly true 
for longer ranged interactions. 

For completeness, we now undertake a brief final anal- 
ysis of a bubble of gas enclosed by a liquid reservoir, 
where Sfi < 0, which is the atipodal system to a drop of 
liquid. The mean-field FMT-DFT approach described in 
Section IV. D is again employed to determine the density 
profiles of the bubble and the intcrfacial properties such 
as the curvature dependence of the surface tension and 
the Tolman length. This is equivalent to the recent DFT 
study reported by Block et al. [73]. As shown in Figure 
9, the curvature dependence of the surface tension for a 
bubble is of a qualitatively different form to that of a liq- 
uid drop, since in this case the dependence is monotonic, 
such that j(R) < 7oo for all R. Taking into account that 
the curvature is now negative R < for a bubble, this 
is consistent with a negative value of the Tolman length 
as for the liquid drop; however, in the case of bubbles 
5(R) remains negative throughout. This general result is 
supported by the findings of recent simulation studies by 
Block et al. [73] . Nevertheless, in other quite recent sim- 
ulation studies employing a mechanical route to analyze 
the data it has been reported that the surface tension 
increases slightly as the radius decreases [55], or that no 
curvature effects for bubbles can be detected [85] . We be- 
lieve that this qualitative discrepancy is again due to the 
inadequacy of a mechanical route to the surface tension. 



VI. SUMMARY AND CONCLUSION 

The purpose of our paper has been to give a compre- 
hensive and up-to-date review of the different approaches 
to the description of structure and interfacial proper- 
ties of microscopic liquid drops and gas bubbles, comple- 
mented with novel mechanical and thermodynamic devel- 
opments, a thorough analysis, and detailed calculations. 
This includes a state-of-the-art description with a non- 
local density functional theory which, in our opinion, is 
the most direct and rigorous way to understanding and 
describing the properties of nanoscale of arbitrary size as 
the approach enables one to represent the marked inho- 
mogeneities of the system. Following a historical intro- 
duction, we started our analysis with a purely mechanical 





FIG. 8. Upper panel: The curvature dependence of the sur- 
face tension of liquid drops obtained with our canonical mean- 
field FMT-DFT for Lennard-Jones particles with different val- 
ues of the cutoff of the potential. The arrows indicate where 
7(J?) = 700 . The values for the planar limit of the surface ten- 
sion are denoted by the dotted lines. Lower panel: The values 
of the planar vapour-liquid surface tension 700 corresponding 
to LJ systems with different cutoffs. The state corresponds 
to a temperature of T* = 0.8 in all cases. 



approach, treating the fluid as a static ensemble of inter- 
acting particles distributed uniformly within the liquid 
phase (of the drop) , while the density of the surrounding 
gas is neglected; the interface is thus perfectly sharp in 
this case (often referred to as the Fowler approximation 
|163j V Such an approach is clearly rather crude but, as 
we show, a simple representation can be developed for the 
vapour-liquid surface tension as the work per unit area 
needed to separate the liquid drop from the rest of the 
fluid. In this type of macroscopic mechanical description 
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FIG. 9. Comparison of the curvature dependence of the sur- 
face tension of a drop and a bubble obtained with our canoni- 
cal mean-field FMT-DFT. The value of the planar limit of the 
surface tension is denoted by the dashed line. The LJ system 
is at a temperature of T* — 1.0 in both cases. 



the surface tension of the spherical interface turns out to 
be proportional to the third moment of the pair potential 
energy; the expression is consistent with that obtained 
in the planar limit by Laplace and Rayleigh in the nine- 
teenth century (cf. Ref. [3 ). In addition we develop a 
novel analytical expression for the Tolman length, as the 
ratio of the fifth and fourth moments of the pair poten- 
tial; the latter measure which characterizes the curvature 
dependence of the interfacial free energy is found to be 
negative from this static mechanical perspective, with a 
magnitude, for typical intermolecular models of simple 
fluids, of about a half of the molecular size. We then 
turned our attention to a purely (macroscopic) thermo- 
dynamic approach that sacrifices the molecular view but 
now incorporates the concept of entropy. Macroscopic 
thermodynamic approaches, as originally introduced by 
Gibbs and then elaborated by Tolman and others, lead 
to a mathematically rigorous description of a liquid drop 
but, as shown in the last part of Section III, can not 
in themselves be used to determine the curvature depen- 
dence of the surface tension neither directly (which would 
require a knowledge of the free-energy density of the en- 
tire inhomogeneous system), nor from the determination 
of the surface of tension. The latter is found to coincide 
with the Gibbs dividing surface so that the surface ten- 
sion takes on its value in the planar limit j(R) — 7oo for 
any R, which would correspond to a Tolman length of 
zero (5 = 0), rendering the Tolman theory inapplicable. 
In order to make progress, the microscopic methods of 
the statistical mechanics are clearly required. 

There are essentially two ways of defining the surface 
tension within a statistical mechanical treatment. One 
approach relies on a mechanical definition of the surface 



tension as the stress transmitted across a strip of unit 
length normal to the interface. In this way, a connec- 
tion between the macroscopic theory of elasticity and the 
components of a microscopic pressure tensor is made to 
determine the surface tension. Beyond the planar limit, 
however, one encounters conceptual difficulties with this 
approach, due to arbitrary nature of the definition of the 
pressure tensor. This lack of uniqueness was already ap- 
preciated by Irving and Kirkwood |Mj and by Harasima 
[98], and later analyzed in a detail for curved interfaces 
by Schofield and Henderson [5] who attributed this ar- 
bitrariness to the fundamental problem of a local def- 
inition of thermodynamic functions depending on two- 
or higher-body interactions (with the exception of the 
chemical potential) . Surprisingly, these warnings are still 
often ignored and a mechanical treatment relying on a 
calculation of the surface tension via the pressure ten- 
sor is frequently adopted. There are two main reasons 
for the ill-advised popularity of a mechanical treatment: 
first of all, any problems related to the non-unique def- 
inition of the microscopic pressure tensor are apparent 
only beyond the planar limit, and, secondly, the use of 
this route is tempting in simulation studies as one re- 
quires a knowledge of the forces between two interacting 
particles in molecular dynamics simulation, and thus the 
components of the pressure tensor can be obtained easily. 

The alternative route is based on the thermody- 
namic definition of the surface tension as the isothermal- 
isochoric change in the free energy per unit area due to 
the deformation of the interface. Within this approach 
one only deals with transformations of the partition func- 
tion so that any problems related with the pressure tensor 
inherent in the mechanical route can be avoided. Never- 
theless, a consideration of first-order changes in surface 
area leads to an expression involving the gradient of the 
potential energy and the correlation function, which for 
a pair-wise interaction gives an identical "stress-strain" 
relation to that obtained with the mechanical route. By 
incorporating an auxiliary external field (the magnitude 
of which is eventually taken to zero) the second-order 
changes in surface area can be analyzed in two ways: 
as a calculation of the magnitude of the external field 
needed to bend the surface; and as the change in the free 
energy accompanying an increase in surface area caused 
by capillary wave density fluctuations. Both methods 
lead to identical formulae involving the one-body den- 
sity distribution function and the direct correlation func- 
tion. These expressions are not exact but sound argu- 
ments provide support for their adequacy up to first or- 
der in curvature and, in particular, its consistency with 
the thermodynamic expression for the Tolman length. 

From a more general point of view, the statistical me- 
chanical expressions relating the macroscopic properties 
of fluids to a given microscopic model can be divided into 
the so-called virial and compressibility approaches. Ex- 
pressions involving the direct correlation function clearly 
correspond to the latter, as commonly implemented in 
a standard statistical mechanical treatment of fluid sys- 
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terns. The use of the term 'virial' can however lead to 
some confusion. It is important to realize that what is 
often referred to as the virial route, is actually only its 
first-order formulation, stemming from the virial theo- 
rem of Clausius |177j . This is the case with the common 
microscopic representations of the mechanical approach 
which amount to a first-order change in free energy per 
unit area implicit in the thermodynamic approach, and 
which have been shown to be valid for homogenous bulk 
systems. On the other hand the virial route is not in 
principle restricted to a first-order representation, even 
though the resulting extension to higher order in the de- 
formation leads to the expressions involving three- and 
higher-body correlation functions. Nevertheless, these 
higher-order terms can be extracted from simulation data 
for distributions of the change in free energy associated 
with the deformation of the interface [72]: the second- 
order (fluctuation) term in the expression for the free- 
energy change turns out to be of the same order of mag- 
nitude as the first-order term for nanoscale drops. This 
clearly emphasizes the fact that the use of a first-order 
virial expression, such as that resulting from a mechan- 
ical treatment, neglects important contributions due to 
fluctuations. This will be the subject of a separate de- 
tailed study [172]. 

All the painful scrutiny and inconsistency of a mechan- 
ical (or virial) treatment can be avoided with the help of 
density functional theory (DFT) , a thermodynamic path 
closely related to that of the compressibility route. In- 
stead of the determination of the direct correlation func- 
tion, however, one simply minimizes the grand poten- 
tial functional to find the equilibrium density profile, and 
then all of the thermodynamic properties of the inhomo- 
geneous system that are required for a Gibbsian thermo- 
dynamic description are available. In our current paper 
we have demonstrated the capability and tractability of 
DFT in providing an unambiguous description of the den- 
sity profile and interfacial properties of liquid drops and 
gas bubbles. Using a non-local mean-field DFT in the 
canonical ensemble together with a consistent Gibbsian 
thermodynamic analysis we come to the following conclu- 
sions for liquid drops: the curvature dependence of the 
vapour-liquid interfacial tension of nanoscopic drops is 
non-monotonic, rising over the value for the planar limit, 
and then decaying slowly to this limit as the radius of the 
drop is further increased; this is consistent with a nega- 
tive Tolman length which we estimate to be about a tenth 
of the molecular diameter; the non-monotonic behaviour 
of the surface tension is reflected in the behaviour of the 
density at the centre of the drop which is seen to cross the 
saturation values of the bulk system at higher tempera- 



tures; our analysis supports the validity of a first-order 
curvature dependence of the surface tension as predicted 
by Tolman relation for drops with microscopic radii down 
to about 10 diameters (below which such a macroscopic 
approach is not expected to be valid); for smaller drops 
it appears that an additional curvature dependence of 
the 1/R 3 form is required in the Tolman treatment of 
the surface tension in order to reproduce the full DFT 
results. The findings for nanoscale bubbles of vapour in 
a bulk liquid are more widely accepted: the curvature 
dependence of the surface tension is monotonic remain- 
ing below the planar limit for all bubble radii; this again 
corresponds to a negative Tolman length which indicates 
that as expected the tension acts on the liquid side of the 
interface. 

Regarding the particular choice of density functional 
(local or non-local, mean-field etc.), it is not yet fully 
clear what impact a given approximation has on the inter- 
facial properties of microscopic drops. It appears, how- 
ever, that to first-order in the curvature of the drop the 
qualitative conclusions are rather insensitive to particu- 
lar form of free energy functional. There now appears 
to be some consensus in the most recent DFT studies 
that the magnitude of the Tolman length is of the order 
of a tenth of the molecular diameter and of a negative 
sign. It is likely, however, that in attempts to go beyond 
a first-order curvature correction, the non-locality of the 
density functional will play a significant role. In studies 
of more complex fluids such as charged particles, poly- 
mer or surfactant solutions, that may exhibit self assem- 
bly one would also expect that a non-local DFT which 
goes further than a mean-field treatment of the attrac- 
tive perturbation term will provide a more appropriate 
description. 
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